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ABSTRACT 


Lifetimes  for  the  1288  KeV  and  the  1440  KeV 


53 

levels  of  Mn  were  measured  using  the  Doppler  shift 

53  53 

attenuation  method  through  the  Cr  (p,  ny)  Mn 


reaction  with  a  proton  energy  of  5.20  MeV.  The 


lifetime  of  the  1288  KeV  level  was  determined  to  be 


p.s.  and  that  of  the  1440  KeV  level  was  found 


p .  s 


The  Doppler  shift  attenuation  method  is  discussed 


in  Chapter  1  with  the  evaluation  of  the  F(t)  function. 

The  mechanism  of  pulse  shape  variation  with  incident 

radiation  in  liquid  organic  scintillators  and  pulse 

shape  distributions  in  Ge  (Li)  detectors  with  their 

timing  effects  are  discussed  in  Chapter  2.  In  Chapter  3 

the  electronics  system  and  experimental  arrangement  are 

explained.  Chapter  4  examines  the  method  of  data  analysis 

53 

and  the  results  of  the  Mn  lifetime  experiment  are  pre¬ 
sented  and  discussed  in  Chapter  5. 
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CHAPTER  1 


THE  DOPPLER  SHIFT  ATTENUATION 
MEASUREMENT  METHOD  (DSAM) 

1.1  The  principle  behind  the  determination  of  nuclear 

lifetimes  by  studying  the  Doppler  shift  or  broadening  of 

gamma  rays  from  the  de-excitation  of  nuclear  levels  is 

based  on  the  rapid  slowing  down  of  moving  atoms  in  a 

solid.  These  swiftly  moving  atoms  have  a  characteristic 

-13 

slowing  down  time  of  the  order  of  10  .  sec.  m  most 

solids  (that  is,  the  time  required  for  the  velocity  of 
the  moving  atom  to  decrease  to  one-half  its  initial  value) . 
The  nature  of  the  function  of  velocity  with  respect  to 
time  is  a  decreasing  exponential,  but  the  explicit  func¬ 
tion  can  only  be  determined  through  a  detailed  examination 
of  the  processes  involved  in  the  stopping  of  atoms  or  ions 
in  solid  materials.  Having  done  this,  however,  one  has  a 
clock  by  which  one  can  compare  the  lifetimes  of  nuclear 

levels.  At  present  the  method  can  be  used  to  determine 

-11  -14 

lifetimes  of  the  order  10  sec.  to  10  sec.  and  is 
limited  by  the  ability  of  the  detection  system  to  resolve 
small  Doppler  shifts. 

The  energy  of  a  gamma  ray  from  the  de-excitation 
of  a  nucleus  moving  at  a  velocity  v  with  respect  to  an 
observer  is  given  by 


3  om  j  e  “  )l  5c  i»r**n  :j,  *i:o  ff  L  n  'ofc  gi 


2 


E  =  E  (1  +  -  cos  0) 
Y  yo  c 


(1) 


or 


AE  =  E  —  cos  0 
y  y°  G 


(2) 


which  is  accurate  to  first  order  in  —  since  the  veloci- 

c 

ties  here  are  less  than  1%  of  c,  the  velocity  of  light. 

E^q  is  the  energy  of  the  unshifted  gamma  ray  and  0  is 
the  instantaneous  direction  of  travel  at  de-excitation 
and  is  measured  from  the  initial  beam  direction. 

Suppose,  now,  the  situation  is  as  depicted  in 
figure  1  where  a  proton  in  the  beam  collides  with  a 
target  nucleus  which  ejects  a  neutron  and  forms  a  recoil 
nucleus  in  an  excited  state.  If  the  recoil  nucleus  decays 
by  emitting  a  gamma  ray  after  a  time  x  (lifetime)  then  it 
will  have  slowed  down  from  some  initial  recoil  velocity 
vq  to  some  velocity  v^  by  encounters  with  atoms  m  the 
target  or  backing  material.  Clearly  v,  =  v-,(x). 


v 


Therefore , 


AE  =  E 

y  y° 


cos  0 


(3) 


or 


AE  =  E 

Y  y°  v 


v,  v 
d  o 


cos  0 
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(4) 


=  E 


v  V  -i 
o  d 


yo 


v 


cos  (0) 


o 


(5) 


V 
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for  the  simple  case  of  a  straight-on  recoil.  E^q  c 
is  the  maximum  possible  Doppler  shift  for  an  infinitely 


3 


Figure  1.  A  proton  (p)  enters  the  target 

material  (in  this  case  also  the  stopping 

53 

material)  and  strikes  a  Cr  nucleus 

forming  a  compound  nucleus  which  decays 

53 

into  a  neutron  (n)  and  a  Mn  recoil 
nucleus  (A)  in  an  excited  state.  The 
recoil  nucleus  then  slows  down  in  the 
stopping  material. 
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Figure  2.  v(t)  is  the  velocity  of  the 
recoil  nucleus  (located  at  the  origin) . 

2  is  the  vector  pointing  in  the  direction 
of  the  gamma  ray  detector.  The  polar 
angles  for  v(t)  and  2  respectively  are 
@2  and  0^,  the  azimuthal  angles  being 
<f>2  and  (j>^.  vq  is  the  initial  velocity 
of  the  recoil  nucleus. 


detector 


FIGURE  2 
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short  nuclear  lifetime  and  one  defines  AE  HE  — . 

yo  yo  c 

The  ratio,  then, 


AE 

AE 


1 

yo 


vd(T) 


V 

o 


F  (t) 


(6) 


is  the  attenuation  resulting  from  the  slowing  down  of 
the  recoil  nucleus,  and  is  equal  to  the  observed  shift 
divided  by  the  maximum  possible  shift. 

In  general,  the  situation  is  more  complicated 
and  the  actual  shift  is  determined  not  only  by  the  speed 
v^(x)  but  by  the  angle  of  recoil.  The  scattered  nuclei 
are  distributed  symmetrically  about  the  initial  beam 
direction  (z  axis)  and, hence,  the  averaged  direction  of 
recoil  is  parallel  to  the  z  axis. 

One  knows  that  a  nucleus  recoiling  in  some  excited 

state  of  lifetime  x  decays  with  a  probability  proportional 

,  -t/x 
to  e  '  . 

It  is  assumed  that  the  recoil  ion  travels  in  an 
initial  direction  denoted  by  the  z  axis  and  that  at  the 
time  of  emission  of  a  gamma  ray  in  a  direction  3  the 
recoil  ion  is  travelling  with  a  velocity  v(t)  making  an 
angle  a  with  3.  3  is  described  by  the  polar  angle  0-^ 

and  the  azimuthal  angle  (j)^  (which  is  zero  here  because 
of  the  choice  of  the  coordinate  system) .  The  direction 
of  v  is  specified  by  0 ^  and  <p^  (figure  2)  . 


t 


The  projection  of  v(t)  on  the  line  of  observa¬ 
tion  is  v(t)  .  3  or  v(t)  cos  a(t) .  Therefore,  the 
Doppler  shift  is  given  by  the  integral  of  this  product 
weighted  by  the  probability  for  decay  of  the  nucleus 
over  all  time  (Ro  69): 

E  OO 

AE  =  ——  /  v(t)  cos  a(t)  exp  (-t/x)  dt  .  (7) 

Y  c  o 


For  a  number  of  such  events,  j,  all 


and  speed  v(t),  and  each  with  angle  (t) , 
Doppler  shift  (over  N  such  events)  becomes 


of  liftime 
the  average 


AE 


Y 


1 

N 


E  N  ~ 

— £  /  v  (t)  cos  a-  exp  (-t/x)  dt 

C  j=l  0  J 


(8) 


Using  the  fact  that  (cos  a)  =  cos  a  and  the 
addition  theorem  of  spherical  harmonics; 

VCOS  “>  =  2lft  J  ,  Yf  (6 l'*l>  Y£(e2'h)  ' 

M=-£ 

one  has  for  the  average  Doppler  shift 


AE 


ao 


Y 


Nc 


N  oo 

£  /  v(t)  ( 

j=l  o 


4  XT 


l 


M*  M 

2H+P  m!_,  yi  (ei'0)  W*’ 

i  i  Xj 


x  exp  (-t/x)  dt 


(9) 
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where  the  values  <f>^  =  0  and  <j>2  =  <j>  have  been  used  (<j> 
is  randomly  distributed  by  reason  of  axial  symmetry) . 

The  probability,  P,  of  an  event  occurring  in  the 
interval  between  <j>j  and  <j>j  +  d<j>  is  proportional  to  d<p, 
or,  if  P  is  normalized  to  unity  over  the  full  range  of 
<f>  (=2tt)  , 


2tt 

P  =  d  <j>/(/  d(J) )  =  d<j>/2  it 
0 


(10) 


Therefore,  for  very  large  N,  corresponding  to  a 
large  number  of  nuclei  decaying  at  the  same  time  giving 
off  gamma  rays  of  energy  E^q  in  their  rest  frames,  one 
has : 


-  =  d  cf)/ 2  tt 


N 


and,  hence,  the  summation  —  £ 

x  2  nr  N  j=l 


(11) 


may  be  replaced  by  the 


integral  /  dd> 
2  7T  /-> 


0 


Expression  (9),  therefore  becomes 


E  00  2  27T  m  m* 

AE__  =  f  v  ( t )  (^4 — )  f  E  jL  ( O 2  (t)  /  $ )  yt  (0t  /0)  dcj) 


Y 


0 


2  £+1  q 


M 


1 


x  exp  (-t/x)  dt 


(12) 


Using  the  well-known  formulae. 


y"(6,$)  =  (-i)-[±|X±  -  p-(cos  0)  exp  (iM<j>) 


M  r  2  £+1  (£-M)'% 


m  j  4jf9  §■ 
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and,  Y”*(6,<j>)  =  (-1)M  If  £M  ( 0 ,  <f, ) 
one  has  for  (12) 


AE 


Y 


Y£ 

c 


oo 

/ 

0 


v  (t) 


x 


( 


2£+l 


) 


2  7T  l 

f  1 
0  M=-£ 


2£+l 
4  7T 


exp  ( iMcJ) )  dc|)  exp 


-M  ,  M 

Pj_  (cos  01)p1(cos  e2(t)) 

(-t/x)  dt  .  (13) 


Clearly,  if  AE  is  to  be  non-zero  then  M  =  0  and  one 

yo 

has 


E 

-  yo 


AE  =  ^  P1  (cos  0,)  f  v(t)  P-,  (cos  0  (t)  )  exp (-t/x)  dt 

y  c  ±  ±  q  J-  z 


E 

-  yo 


oo 


cos  0 /  v(t)  cos  09  (t)  exp  (-t/T)  dt  .  (14) 

1  0 


Of  course,  initially  the  recoil  ions  are  moving 


-»■ 


with  velocity  vq  (along  the  z  axis)  and  so  the  maximum 
Doppler  shift  for  these  ions  are  given  by 


E  «> 

AE  =  — ^  /  v  .3  exp  (-t/x)  dt 

yO  r*  r\ 


c  0  ° 


U5) 


using  (7)  with  v(t)  =  v  ,  a  constant.  Or 
* 

E 

AE  =  — ^  /  v  cos  0,  exp  (-t/x)  dt  (16) 

yo  c  q  o  1 

where  the  expression  vq.3  =  vq  cos  0^  has  been  used. 
From  the  previous  definition  of  F(x)  one  has: 


.  ,* 
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F  (x) 


00 

f  v(t)  cos  <j)(t)  exp  (-t/x)  dt 
0 

00 

/  VQ  exp  (-t/x)  dt 


(17) 


where  (t)  =  4)  (t)  .  According  to  convention  one  defines 

cos  cp  (t)  to  be  cos  <j>  (t)  which  projects  the  average 
velocity  v(t)  onto  the  z  axis  which  is  the  initial 
direction  of  travel. 

1.2  In  order  to  determine  the  velocity  function,  v(t), 

dE 

it  is  necessary  to  know  the  stopping  power  ^  of  the 
slowing-down  material.  E,  of  course,  is  the  energy  of 

f. 

the  recoiling  ion  and  x  is  its  distance  traversed. 

From  simple  mechanics  one  has 

E  =  P2/2m  (18) 


where  P  is  the  momentum  of  the  ion  and  M  its  mass. 
Therefore , 


dE  (v) 
dx 


I  dP 
M  dx 


M 


dv 

dx 


v 


=  M 


dv  dx 
dx  *  dt 


=  M 


dv  (t) 
dt 


(19) 


Hence , 

M  dv 
dE  (v) 
dx 


dt 


(20) 


Jc.  .  ei3M  ci  .  Vo  3  .UUsanm  m»6-eftb»Je 
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can  be  used  to  derive  v  (t)  .  — -  is  usually  given  in 

2 

units  of  KeV.cm  /mg  and  since ,  dimensionally, 

-  ]  2-1  - ^ 

MeV.cm  ‘  =  (KeV.cm  .mg  )  x  (g.cm  ) 

(the  last  factor  on  the  right  is  clearly  the  density 
of  the  material,  p) 


dE 
dx ' 


(MeV .cm  ^ )  = 


d,E 

dx 


(KeV.cm  .mg 


(21) 


and  x  =  x 1  p 


(22) 


in  transferring  from  one  system  of  units  to  another. 
dE 


dx 


can  be  broken  down  into  the  sum  of  two  distinct 


processes:  electronic  and  nuclear  stopping.  That  is, 


dE  _  .dE.  ,dE. 

dx  dx  e  ‘  dx  n 


(23) 


where  the  subscripts  e  and  n  refer  to  the  electronic 
and  nuclear  processes  respectively. 

Electronic  stopping  ,  the  mechanism  of  which  is 
essentially  the  collision  between  electrons  from  the 
surrounding  medium  and  the  recoiling  ion,  is  a  con¬ 
tinuous  process  since  each  individual  encounter  with  the 
electron  cloud  involves  an  energy  loss  of  the  order  10  z^ 
electron  volts  (No  63)  ,  where  is  the  atomic  number  of 
the  stopping  atom,  which  is  much  smaller  than  the  recoil 


11 


energy  of  the  nucleus  which  is  of  the  order  of  0.1  MeV. 
Hence,  this  process  does  not  involve  changes  in  the 
recoil  direction.  Nuclear  encounters,  which  take  place 
at  much  lower  energies  than  the  initial  recoil  energy 

vc 

(137  <<  1,  where  v  is  the  velocity  of  the  ion  when 

c 

nuclear  stopping  becomes  larger  than  electronic  stopping) 

(No  63),  involve  large  direction  changes  which  are 

manifested  by  the  cos  <}>  (t)  term  in  the  expression  (17) 

for  AE  . 

Y 

1.3  Electronic  Stopping 

One  can  gain  insight  into  the  behavior  of  the 
electronic  stopping  process  involving  the  encounter 
between  a  moving  ion  and  the  electron  cloud  surrounding 
the  atoms  of  the  stopping  material  by  considering  this 
cloud  to  be  a  degenerate  electron  gas.  That  is,  the 
electrons  (Fermions)  are  assumed  to  fill  all  available 
states  up  to  a  certain  energy  (yF) ,  with  all  states  of 
energy  greater  than  y^  being  completely  empty.  The 
velocity  of  the  moving  ion  (V)  is  assumed  to  be  consi¬ 
derably  less  than  the  maximum  allowable  electron  velocity, 

v^  (V  «  v„)  .  Hence,  for  an  electron-ion  collision  the 
F  F 

velocity  change  suffered  by  the  electron  will  be  of  the 
order  of  V  (the  maximum  possible  change  will  of  course 
be  2V) .  Clearly,  then,  collisions  for  which  the  final 


■ 


■ 
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electron  velocity  lies  inside  the  occupied  zone  of  the 

electron  momentum  space  will  be  forbidden  by  the  Pauli 

principle.  The  only  collisions  allowed  are  those  with 

electrons  having  a  velocity  differing  from  v^  by  a  value 

r 

of  the  order  of  V.  In  momentum  space  the  electrons 
occupy  a  volume  of  dimension  P  =  M  v„,  where  M  is  the 

r  r 

mass  of  an  electron.  The  electrons  capable  of  collision 

occupy  a  thin  shell  between  P^  and  Pp  +  dP^  so  that  the 

2 

volume  formed  by  this  shell  will  be  approximately  P^dPp. 

Since  the  minimum  volume  in  momentum  space  that  can  be 

3 

occupied  by  a  single  electron  is  fi  ,  the  number  of  electrons 
engaging  in  collisions  will  be  given  by 


n 


PF  dPF 


(24) 


But  by  previous  argument  dP^  ~  mV  so  that  one  has 


n 


(25) 


The  interaction  between  the  ion  and  the  electron 
is  basically  Coulombic  so  the  cross-section  for  collision, 
a,  will  be  proportional  to  the  Rutherford  scattering 
cross-section.  Therefore 


<^2> 

M  v„ 


a 


2 


(26) 
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where  z  is  the  atomic  number  of  the  ion. 

In  this  type  of  collision  not  forbidden  by  the 
Pauli  principle  the  energy  transferred  during  a  colli¬ 
sion,  W,  will  be  given  by 

W  s  j  M  [  (vp  +  V)  -  vp] 

~  M  vp  V  (27) 

to  first  order  in  V. 

The  total  number  of  collisions  taking  place  per 
unit  time  will  be  n  avp,  and  hence,  the  energy  transfer 
per  unit  time  will  be  W  n a .  Thus  one  has  for  the 
energy  loss  per  unit  time  for  the  ion 


or, 


Then, 


dW 

dt  =  W  n  a  VF 


dW  dX 
dX  dt 


M2e4 

-ft3 


V' 


dW 

dX 


V  = 


m2  4  0 

Me  v2 

ftJ 


(28) 

(29) 

(30) 


So,  the  energy  loss  per  unit  distance  travelled  by  the 
ion  is  given  by 


dW 

dX 


cc  V 


(31) 


For  this  very  simple  case,  then,  it  has  been 
shown  that  the  stopping  power  due  to  electronic  slowing 


. 
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down  is  indeed  proportional  to  the  velocity  of  the  ion. 

In  actual  fact,  the  process  is  very  complicated  and  one 
is  required  to  determine  the  probability  that  an  electron 
will  be  captured  and  then  lost  by  the  moving  ion.  When 
the  ion's  velocity  is  much  greater  than  the  orbital 
velocities  of  the  electrons  its  electrons  are  almost 
entirely  stripped  away  upon  entering  slowing-down  material. 
Initially,  the  ion  loses  energy  by  collisions  with  the 
electrons  of  the  material  which  will  not  be  retained  by 
the  ion  because  although  there  is  a  small,  finite  proba¬ 
bility  of  capture  of  an  electron  during  one  collision 
there  is  a  very  much  larger  probability  of  its  loss  during 
the  next  collision.  As  the  velocity  of  the  ion  decreases, 
however,  to  velocities  of  the  order  of  the  electron  orbital 
velocity,  the  probability  of  electron  loss  decreases  and 
that  of  electron  capture  increases  until,  finally,  at 
ion  velocities  below  the  orbital  velocity  the  loss  pro¬ 
bability  approaches  zero  with  a  further  increase  in  the 
capture  probability.  At  this  point  an  electron  is  cap¬ 
tured  and  held  by  the  ion  in  one  of  the  lowest  lying 
atomic  shells.  To  capture  a  second  electron  in  a  higher 
shell  (with  a  correspondingly  lower  orbital  velocity  than 
for  the  first  shell)  the  ion  must  slow  down  to  the  next 
highest  orbital  velocity  by  the  scattering  process 
discussed  above.  Once  the  second  electron  is  captured 
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a  third  may  be  captured  (it  will  have  a  lower  orbital 
velocity  than  either  first  or  the  second  electron  since 
it  will  be  captured  into  a  "higher"  level  than  either 
of  these)  in  a  similar  manner  and  so  on  for  all  the 
electrons  (z)  capable  of  being  held  by  the  ion,  the 
only  difference  for  each  capture  being  the  necessity 
of  the  ion  to  slow  to  the  successively  lower  orbital 
velocities  of  these  electrons  within  the  ion. 

At  the  point  when  the  ion  has  slowed  down  to 
velocities  less  than  the  orbital  velocity  of  the  least 
tightly  bound  electron,  the  ion  is  mostly  a  neutral 
atom,  and  further  energy  loss  proceeds  predominantly 
by  the  elastic  collisions  of  the  screened  electric 
fields  of  ion  and  stopping  nuclei.  The  lower  the 

velocity  of  the  ion  the  more  the  latter  process  dominates, 

> 

and  the  less  is  the  energy  loss  to  the  orbital  electrons. 
The  final  result  is  of  course  the  stopping  of  the  moving 
ion;  that  is,  the  thermalization  of  the  ionic  velocity. 

Note  that  in  the  simple  case  considered  the 
velocity  of  electrons  above  the  conduction  band,  vp, 

o 

is  approximately  8  x  10  cm/sec  using  Morse's  expression 
(Mo  64)  with  p  =  10  and  that  the  velocity  of  electrons 

g 

in  the  Bohr  atom  is  2.18  x  10  cm/sec.  So  the  condition 
V  <<  v^  would  indicate  that  at  these  velocities  the  ion 
would  already  have  captured  electrons  into  the  lower 
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atomic  shells  before  the  simple  case  could  even  be 
considered.  All  this  picture  has  shown  is  the  appro¬ 
ximate  linear  dependence  of  the  stopping  power  on 
velocity,  a  fact  which  has  experimentally  been  verified 
(Po  61) ,  (Re  55) ,  and  (No  60) . 

Lindhard  and  Scharff  (Li  61) ,  (Li  63)  have  used 
the  Thomas-Fermi  statistical  model  of  the  atom  to 
obtain  a  theoretical  expression  for  the  stopping  power: 


dE 

dX 


K 


E 

M 


KeV  cm  /mg 


(32) 


where  E (v)  is  the  energy  of  the  moving  ion,  is  its 
mass  and 


73.0  z 


1/6 


K  = 


M2(z2/3  + 


2/3  v  3/2 

z2  > 


for  v 


$ 


fi 


2/3 


(33) 


z ^  and  z ^  are  the  atomic  numbers  of  the  moving  ion  and 

the  stopping  material  atom  respectively  and  M2  is  the 

mass  of  the  stopping  atom. 

Firsov  (Fi  58) ,  (Fi  58) ,  and  (Fi  59)  proposed 

an  alternate  method  of  deriving  a  theoretical  approxi- 
•« 

mation  for  the  stopping  power.  He  considered  the 
regions  of  the  interaction  potentials  of  the  two  colli¬ 
ding  ions  to  be  defined  by  a  surface  S,  so  that  the 
average  effect  of  electrons  crossing  S  from  one  region 
to  another  is  that  they  assume  the  momentum  of  the  atom 
in  the  latter  region. 


* 
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The  energy  loss  of  the  electrons  was  given  by 
Firsov  as: 

**■  • 

£  =  M  f f  —•  dS  (34) 

S 

where  n  and  v  are  the  number  density  and  the  velocity 
of  the  electrons  respectively  before  crossing,  &  is  the 
radial  vector  between  the  two  nuclei,  and  d§  is  a  surface 
element  of  S.  v  is  then  written  in  terms  of  n,  n  being 
written  in  terms  of  the  potential  <j>  by  considering  the 
Thomas-Fermi  model  of  the  atom.  Therefore, 

M2  2  9 

€  =  "9e-~  ff  (f)  dS  R  d  R  .  (35) 

4tt  fi  S 


as : 


The  interaction  potential  was  given  by  Firsov 


<i>  (r) 


X[^(V  Z2]  !] 


(36) 


where  x(x)  is  screening  function  in  the  Thomas- 

2/3  2/3 

Fermi  potential  and  where  a  =  aQ  0.8853  (z^  +  z ^  ) 

is  the  screening  length,  aQ  being  the  Bohr  radius.  The 
function  ip  =  (z^3  +  was  assumed  to  be  true 

over  most  of  the  range  of  the  potential. 

From  these  relations  Firsov  obtained 


dE 

dX 


19.71  ( z  x  + 


Z2) 


M, 


(JL) 

VM.  ' 


(37) 
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for  the  stopping  power  which,  like  Lindhard's  formula 
and  the  very  crude  picture,  predicts  a  linear  velocity 
dependence . 

These  relations  are  only  approximate,  however, 
and  represent  a  "smoothing"  over  all  electron  orbits. 
Using  ion-gas  (Mo  62) ,  (Af  65) ,  (Ev  65) ,  (Ke  65)  and 
ion-surface  studies  (Ra  64) ,  (Da  65) ,  (Pi  64) ,  (Ro  65) 
it  has  been  shown  that  the  electron  excitation  losses 
occur  in  a  step-wise  fashion,  the  electrons  being  removed 
with  discrete  energies  as  a  result  of  overlapping  of 
successive  energy  levels.  Indeed,  the  stopping  power 
has  been  shown  to  be 


dE 


dx 


where  0.8  <  S  <1.4  (Ro  69),  and  c  is  a  constant  of 
proportionality . 

1.4  Nuclear  Stopping 


In  this  stopping  process  momentum  and  kinetic 


energy  are  transferred  to  the  translationary  motion  of 
the  target  atom  as  if  it  were  a  hard  sphere  (elastic 
collisions).  Thus,  the  problem  is  much  simpler  than 
the  electronic  case  where  energy  was  transferred  to  the 
individual  electrons  of  the  target  atom  producing  atomic 
excitation  and  ionization  (inelastic  collisions) .  The 


. 
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scattering  field,  when  screening  is  considered,  leads 
to  a  potential  of  the  incident  particle  of  the  form 
(Bo  48) 

z  e  -r/a 

<P  (r)  =  -  e  (39) 

where  r  is  the  distance  from  the  center  of  the  atom, 

z^e  and  z^e  are  the  charges  of  the  incident  and  target 

particles  respectively,  and  a  is  the  screening  length 

with  a  =  a  S  \  a  being  the  Bohr  radius  and  S  is  the 
o  o  ^ 

1/3 

screening  factor.  For  a  single  atom  of  charge  ze,  S=z  ' 
is  an  approximate  estimate  of  the  screening  factor.  In 
the  case  of  a  collision  involving  two  atoms  the  two 
screening  factors  may  be  added  in  quadrature  to  obtain 
an  approximate  resultant  factor  S  given  by 

S  =  sTl  +  S2  (40) 

where  S-^  =  z^^  and  =  z\^  are  the  individual'  screen¬ 
ing  factors.  S  roughly  accounts  for  the  interaction  of 
the  two  charge  distributions  corresponding  to  potentials 
of  the  type  given  by  (39) . 

2  2 

Define  the  quantity  E,  =  b/a,  where  b  s2|z^z^e  |/(Mqv  ) 
is  the  collision  parameter  (b  would  be  the  minimum  distance 
of  approach  for  two  repelling  particles),  Mq=M1M2/ (M^+  M2) 
is  the  reduced  mass  of  the  system  of  two  colliding  parti¬ 
cles  of  masses  and  M2 .  The  electron  orbital  speed  of 


. 
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2 

the  hydrogen  atom  is  vq  =  e  /ft.  Let  A  be  the  wavelength 
of  the  particle  of  mass  Mq  and  velocity  v.  Then, 


A 

and  if 

X 


ft 


M  V 
o 


b/A  , 


therefore , 

X  = 


fl  v 


2TZ2 


V 

o 

V 


(41) 


(42) 


and 

5  =  2  z  z2  S  /2)  (43) 

o 

where  y  is  the  mass  of  the  electron. 

g 

Since  vq=  2.18  x  10  cm/sec  and  since  the  maximum 

g 

velocities  that  will  be  considered  are  v  :$  10  cm/sec 
and  both  particles  are  assumed  to  carry  considerable 
charge  (z^,  z ^  *  25)  one  always  has  the  conditions  that 
X  >>  1  and  £  <<  1.  These  conditions  will  remain  over  a 

considerable  range  of  variation  of  the  assumptions  since 

-4 

the  ratio  y/M  will  be  of  the  order  of  10 

o 

Consider  now  a  particle  passing  through  a  layer 
of  matter  Ax  in  thickness  with  N  atoms  per  unit  volume 
(Bo  48) .  Further,  it  is  assumed  that  Ax  is  small  so 
that  the  initial  velocity  of  the  particle  is  not 
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appreciably  changed  upon  its  passage  through  the  thick¬ 
ness.  The  number  of  collisions  is  then  given  by: 

u  =  N  Ax  a  (44) 

where  a  is  the  collision  cross-section. 

This  frequency  of  collisions  will  fluctuate 
statistically  and  the  probability  of  having  n  collisions 
in  an  interval  of  the  momentum  transfer  T  (from  one 
particle  to  the  other)  will  take  on  a  Poisson  distri¬ 
bution  : 

P  (n)  =  ^  e~“  (45) 

II* 


where  the  average  number  of  collisions  in  the  interval 
is  w. 

In  the  i^*1  interval  of  the  total  collision 
process  (formed  by  the  successive  collisions  between 
the  incident  particle  and  the  particles  in  the  thickness 
Ax)  with  kinetic  energy  transfer  T^  and  number  of  colli¬ 
sions  n^  (the  average  number  being  or  )  the  energy  loss 
is  given  by 


n . 

l 


T  . 

l 


and,  over  all  intervals  one  has  that 


(46) 


AE 


l 


(47) 
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Replacing  with  its  mean  value  gives  an  average 
energy  loss 


AE  =  £  T.  co  . 

•  -L  JL 
i 


(48) 


If  these  intervals  are  made  infinitesimally  small 
the  summation  becomes  an  integral  so  that 


AE  =  /  T  dco  =  N  Ax  /  T  da 


(49) 


Ax  being  a  constant  element  of  the  total  range  R. 

da  may  be  evaluated  very  easily  by  considering  the 

simple  case  of  a  Rutherford  type  of  collision: 

2  2 


da 


/  1  2  .  dft 

o  m  2  .40 

2  M  v  sm  tt 
o  2 


(50) 


where  d£3  =  2 tt  sin  0  d0  is  the  solid  angle,  0  is  the 
deflection  suffered  by  the  incident  particle  in  the  center 
of  mass  system  (figure  3a) ,  Mq  is  the  reduced  mass  defined 
above,  and  v  is  the  initial  velocity  of  the  incident 
particle . 

From  the  kinematics  of  the  collision  the  kinetic 
energy  transferred  during  a  single  collision  is  given  by 


T 


.  2 
sm 


0 

2 


where 


4  Mx  M2 

(Mj+  M2)2 


E 


(51) 


Figure  3a.  The  collision  between  particle  1 
moving  from  the  left  with  initial  velocity  v 
(in  lab  system  with  particle  2  initially  at 
rest)  and  particle  2  moving  from  the  right 
with  the  velocity  of  the  center  of  mass  v 

c 

is  shown.  Particles  1  and  2  are  scattered 
through  an  angle  0  in  the  center  of  mass 
system  and  particle  2  is  scattered  through 
angle  0 ^  in  lab  system  acquiring  energy  T. 
Equation  51  follows  immediately  from  this 
diagram. 

Figure  3b.  Hyperbolic  orbits  of  particles 
1  and  2  in  the  center  of  mass  system  are 
shown  with  the  asymptotic  angles  a.  F  is 
the  focus  of  orbit  1  and  r  is  the  vector 
joining  the  focus  with  particle  1  and  makes 
an  angle  0'  with  the  vertical  symmetry  axis. 
The  impact  parameter  P  is  shown. 
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FIG  URE 


2  4 


•  1  2 
is  the  maximum  kinetic  energy  transfer  and  E  =~  v 

Therefore  substituting  (51)  into  (50)  gives 


da  = 


0  2  2  4 

2 it  z^  z 2  e 

M2  v2 


dT 

n,2 


=  B 


n 


dT 

T2 


(52) 


and  hence , 


AE  = 


T 

N  Ax  B  / 
n  T 


M 


dT 

T 


T 

=  N  Ax  Bn  log  T 


M 


(53) 


In  order  to  evaluate  T  ,  the  minimum  kinetic 

a' 

energy  transfer,  consider  the  scattering  case  depicted 
in  figure  3b.  From  the  diagram  one  can  see  that 


2a  +  0  =  7T 


or 


0  _  7T 

2  2  a 


(54) 


and,  therefore 


tan  f 


=  cot  a 


Now  from  the  solution  of  the  orbital  equation 


(Sy  64) 
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M  K 


o 


(55) 


one  has 


u  =  —  =  B  +  A  cos  (0  -  0  ) 
r  o 


(56) 


2 


where  K  =  z^z^e  ,  L  is  the  orbital  angular  momentum. 


and  B  =  -M  K/L2.  Define 


o 


e  =  A/ | B | 

which  is  the  eccentricity  of  the  orbit,  and 

B 


one  finds,  then 

a  (e2-  1) 

zz  - - - - - 

e  cos  0-1 

so  that  r  ->  °°  (0  +  a)  for  cos  a  =  1/e  or. 


cot  a  =  (e2  -  1)  2 


(58) 


which,  by  evaluating  A  from  the  initial  conditions  of 
the  orbit  and  thus  finding  that 


2  E  L 


2 


e 


(1  + 


2 


(59) 


M  K 


o 


. 
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becomes 


tan 


0 

2 


cot  a 


M  K2  % 
—2 _ ) 

2 

2  E  L 


(60) 


since  it  is  assumed  that  momentum  is  conserved  L  =  M 
giving  the  result  (P  is  defined  in  figure  3b) 


o 


v  P, 


tan  | 


(61) 


one  should  observe  that  for  p  small  compared  to  a,  colli¬ 
sions  will  occur  only  in  the  unscreened  part  of  the 

field.  Therefore,  the  minimum  angle  0  will  be  given 

a 

by 


tan 


A 

2a 


(62) 


It  is  permissible  to  accept  only  first  order  terms  in 
ea 

-y  since  the  condition  E,  <<  1  holds. 

However  since 

2  0a  9a  2 

Ta  =  TM  Sln  T  ~  TM  i2-)  (63) 

then , 


£  ■  <V2  -  <+>'2 


2  E  M2  hi 


#  ,  4  ,  2/3  2/3  x  ^ 

(Mi+M2)ziz2e  y(Zj_  +  z2  ) 


(64) 
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Therefore  one  has  for  the  nuclear  stopping  power  upon 
combining  (53)  and  (64) 


=  N  B '  log 
Ax  n  ^ 


2  E  M2  ft' 


A  4  ,  2/3  2/3, k 

_(M1+M2)  z1z2e  u(z1/  +  z27  )2 


(65) 


where 


B  = 
n 


o  2  2  4 

2 tt  z ^  z2  e 

“T1  ‘ 


(66) 


Lindhard  and  Scharff  have  introduced  the  dimension¬ 
less  variables 


4 tt  N  M0  a2  M, 

- ± x  (67) 

(Mx+  M2)Z 

a  M0 

- * -  E  (68) 

z-^  z2  e  (M^+  M2) 


and 


e  = 


with 


a 


(z 


0.8853 

?/3+  Z2/3 


2 

e 


(69) 


Using  these  variables  (65)  becomes 

P-  =  p  log  (2.259  e)  .  (70) 

dp  2e 

According  to  Robertson  (Ro  69)  in  his  evaluation 
of  f(t),  the  nuclear  stopping  formula  due  to  Lindhard 
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(72)  can  be  approximated  by  the  following  functions: 


for  e  <  14  , 

dje  =  1  i-  ( £  -  14)  log  (2.259  e) 
dp  4  2  £ 

+  (18  -  £ )  (1.97  £^  exp  (-1.82  e^) 

+  0.0225  £  2 }  ,  for  14  ^  ^  18 

=  —  log  (2.259  e)  for  £  >  18 

so,  from  the  last  relation  one  sees  that  both  the  appro¬ 
ximation  and  the  Lindhard  stopping  formula  join 
smoothly  with  the  Rutherford  stopping  in  the  high  energy 
region.  The  above  functions  agree  with  (72)  to  within 
20%  over  the  range  of  £  considered. 

1.5  F(t)  Calculation  of  Blaugrund  (B1  66) 

Upon  examining  the  behavior  of  the  recoil  ions 
in  a  stopping  material  it  has  been  seen  that  the  two 
main  processes  of  energy  loss  are  electronic  and  atomic 
collisions,  and  that  atomic  collisions  can  be  further 
divided  into  kinetic  energy  losses  and  atomic  scattering. 
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The  previous  discussions  have  provided  enough 
information  to  determine  v(t),  the  speed  of  the  recoil 
ion  in  the  stopping  material,  but  in  order  to  evaluate 
(17)  it  is  necessary  to  include  the  effects  of  scatter¬ 
ing  which  results  in  the  cos  cj>  (t)  term  in  the  expression 
for  F  (t)  . 

Lindhard  and  Scharff  (Li  63)  derive  a  screened 
potential  (^)  ( r )  by  an  extrapolated  perturbation  method 
for  classical  scattering,  and  from  this  obtain  an  expres¬ 
sion  for  the  atomic  scattering  cross-section: 


da 


tt  a 


dt 


2t 


3/2 


f  (t2) 


(71) 


where  t^  =  *  ® 


e  sin  2  ,  s  being  the  dimensionless  quantity 
given  by  (68),  and  as  before  sin  =  T/T  .  f(t2)  i-s 

called  the  universal  differential  scattering  cross-section 
for  elastic  nuclear  collisions  and  the  behavior  of  this 
function  is  shown  in  figure  2  of  Lindhard' s  paper  (Li  63) 
and  is  derived  using  a  Thomas-Fermi  potential. 

It  follows  from  the  above  definitions  that  t/e  is 
the  dimensionless  kinetic  energy  and  using  the  expression 
da/ra  as  the  dimensionless  cross-section  (49)  becomes 


dp '  n 


|  da  (tb 

Tra  0 


/  f  (x )  dx 
0 


(72) 
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i, 

using  (71),  where  x  =  t2. 

From  before  for  the  low  ion  velocity  range 
e^  2 /3 

(v  <  z^  )  the  electronic  stopping  power  can  be 

written  in  dimensionless  units  as 


(~) 

vdp'e 


=  k 


(73) 


where 


k  =  z 


1/6 


0.0793  z' 


(A. 


+  a2) 


3/2 


,  2/3 ,  2/3, 3/4  3/2 

(Z1  +  Z2  }  A1  A2 


Two  new  parameters  are  introduced: 

v  H  —  is  the  dimensionless  velocity  of  the  ion, 

A  c 
e 

and  771  is  defined  by  the  relation  e  =  \1YLv  ,  so 
that 


1.63  x  10 


A1  A2 


,  2/3 
Z2(Z1  + 


z2/3)!S<A1+  A2> 


(74) 


One  can  also  define  a  dimensionless  time  variable, 
0,  such  that: 


(75) 


Thus,  one  has  6  =  t/T  where  t  is  time  in  seconds  and 
the  coefficient  T  is  given  by 
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T  = 


fi  (Aj+  A2)2 
e2  4tra2  N  A-^  A 


(76) 


Using  these  expressions  the  integral  (20)  becomes 


e  =  (|  Trl )  V 


de 


^  de 


(77) 


dp 


where  de/dp  =  (de/dp.)  n  +  (de/dp)  and  eQ  is  the  energy 
of  the  moving  ion  at  t=0 .  This  integration  may  be 
performed  numerically  to  give  v  as  a  function  of  0. 


Having  defined  these  quantities,  cos  <j> ,  the  average 
projection  of  the  velocity  onto  the  initial  beam  direc¬ 
tion,  may  now  be  evaluated. 

Lewis  (Le  50)  has  given  the  angular  distribution 
for  an  initially  collimated  beam  of  ions  after  passing 
through  a  distance,  x,  of  a  stopping  material: 


00 


X 


F  (<j) , x)  =  j—  Z  (2&+1)  P  (cos  <J>)  exp  (-/  K  dx) 
qiJ  1=0  36  0  9 


where 


(78) 

(79) 


K_  .  =  N  /  [1  -  P0  (cos  <J))  ]  da  (4>)  , 

Z(p  x, 

■% 

N  being  the  number  density  of  stopping  atoms  and 
da(4>)  =  a  (4))  dft^  is  the  differential  cross-section 
as  a  function  of  the  lab  angle,  (p  .  P^  (cos  4>)  are  the 
unnormalized  Legendre  polynomials.  (For  an  outline  of 
the  derivation  of  (78)  see  appendix  1) . 


■ 
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From  (78)  one  gets  immediately , 


P  (cos  d>) 
n 


/  F  (cp  ,x)  Pn(cos  <j>)  dj  d0 


25  +  1  ^Tr  1  p 

~T—  f  /  P  (cos  (p)  d  (cos  (j) )  de 
77  0  -1 


x 


x  exp (-/  K  dx) 


0 


n<j> 


=  exp  (- 


x 
f  K 
0 


n<f> 


dx) 


(80) 


or,  for  n  =  1 


cos  $  =  exp 


x 

(-  /  K 
0 


14> 


dx) 


(81) 


One  also  has  the  well-known  expression  relating 
the  lab  angles  to  the  center  of  mass  angles  (J>  and  0 
respectively: 

tan  4>  =  — Sln  9  ■  (82) 

r  +  cos  0 


where  r  =  A^/Ap,  A ^  being  the  mass  of  the  recoil  ion 
and  Ap  that  of  the  stopping  ion. 

Solving  (82)  for  cos  $  one  finds 

r  +  cos  0 

(r  +  2r  cos  0  +  1)  2 


cos  (p  = 


(83) 


■ 
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which,  when  expanded  in  terms  of  Legendre  polynomials 
(see  appendix  2)  becomes 


00 

'ppcos  0)  +  £  (-l)n  2^TT(Pn+i(cos  0)-Pn_1(cos  0)]rn 


COS  (j)  = 


for  r  <  1 


(84a) 


00 

1  +  Z  (_1)n  2HTT[Pn+l(cos  0)-p„_!  (cos  0)  ]  r~  <n+1) 


n=l 


n-1 


for  r  >  1 


(84b) 


In  expression  (80) 


K£<£  =  N  *  “  p £  (cos  4>)  ]  da  (<J>) 


Putting  the  expression  for  cos  <j>  into  k-^  one  gets 


Kl0  =  N  /  [1  -  cos  <J)  (0)  3  da  (0) 


(85) 


oo 


=  N  /  [1  -  ppcos  0)]  +  £  (-1) n  rn 

n=l 


x  [-Pn+1(cos  0)  +  pn_1(cos  0) ] da (0) ,  (86) 


using  (84a)  and  using  da(0)  instead  of  da((J>)  in  (79). 

The  first  square-bracketed  term  in  (86)  is 
immediately  seen  to  be  k^0,  where 


k10  =  N  f  t1  “  px(cos  0)]  da (0) 


(87) 
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The  second  bracketed  term  may  be  rewritten  in 


the  form. 


[1  -  pn+]_(cos  0)1  “  [1  “  pn_i(GOS  Q))'  (88) 


so  that  now  these  two  terms  taken  into  the  integration 
become 


oo 


Z  (-1) n  rn  (k  .  -  k  ln  )  . 

1  2n+l  n+10  n-10 

n=l 


(89) 


Then, 


X  X 

- r  /  /*  i  x  n  n+1  n 

cos  <j>  =  exp  (-/  k  dx).exp(-/  (-1)  9  r 

0  u  0 


X  [kn+10  kn-10 


(90) 


OO 


p, (cos  0)  n 
1  n=l 


Pn+1(cos  0) 


pn_i (cos  0) 


,  -i  \  n  ,  n+1 .  n 

(~1}  (^+r)r 


r  <  1 


(91) 


Similarly  for  r  >  1 


OO 


cos  <j>  =  n 
n=l 


Pn+l(COS  8) 

P  '  (cos  0) 

n-1 


(-1) 


n+1  n  1 


2n+l  n+1 
r 
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using  the  properties  of  the  exponential  function  and 
the  relation 


x 


P  (cos  0)  =  exp(-  /  k  ^  dx)  . 
nx  '  ^  q  nQ 


In  order  to  calculate  the  P  (cos  0)  terms  it 

. ^ 


will  be  necessary  to  calculate  cos  Q,  which  is  given 


by  (81)  and  (79)  upon  replacing  (j)  with  0  and  setting 
£  =  1.  The  integrand  of  (79)  can  be  written  as 


/ 1  0  .  2  0 
(1  -  cos  0)  -2  sin  2 


(92) 


From  before  we  have 


j  2  dt  £ 

da  (0)  =  ra  -  -  * 


2t 


3/2 


f  (t2) 


where  t"2  =  e  sin  ~  . 


O  r]  f  4- 

da  (0)  =  2  •  ’  ~2  0  f(t2) 

e  sin  2 


=  Tra 


d^ 


2  .20 
e  sm  2 


f  (y)  /  y  =  t' 


(93) 


Therefore 


[1  -  cos  0]  da(0)  = 


2  7T  3.  r  /  \  j 

— 2 —  f  (y)  dy 


•  • 


giving 


K 


10 


2irNa  ,d£* 

£  ldp;n 


(94) 


‘ 
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by  using  (79)  and  since 


(§£) 

dp'  n 


f  (y)  dy 


(72) 


which  is,  in  fact,  an  average  of  f (y)  and  thus  is 
assumed  to  be  a  constant. 

Examining  the  integration  over  x  in  (81)  recall 
the  previous  definition  of  p;  so  one  may  write 


x  = 


47TNa‘ 


(A1+  V 

A1  A2 


P  = 


(r  +  1) 


47rNa‘ 


(95) 


Thus  dx 


“  dp  =  de/(f|)  • 


X 

f  K 
0 


10 


d  x  =  / 


°  O  M  2 

27TNa 


Hence , 
(de/dp) 


n 


(r+1) 


de/dp 


4iTNa' 


de 


1  (r+1) 2  ,°  (d£/dp ) n  _  1  (r+1) 2 

2  r  £  s  (de/dp)  E  ~  2  r 


(96) 


where 

£o  (de/dp)  ^ 
1  5  1  e  (de/dp) 


Therefore , 

2 

cos  0  =  exp  [-  J  ^  -  I]  •  (97) 

For  the  case  r  =  1  one  has  <j>  =  j  0  and  cos  <f> 
may,  in  principle,  be  calculated  directly  without  any 
approximations.  Using  the  approximation. 
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P^(cos  0)  ~  (cos  0) 


l 


(98) 


one  has  for  r  ^  1 


/  -i  \  n  2  (n+1)  n  .. 

(_1)  -5H+r- r  '  r  < 1 


cos  cj>  = 


cos  0  n  [cos  0] 
n=l 


(99a) 


oo 


n  [cos  0] 
n=l 


,  i v  n+1  2n  1  n 

(_1)  2n+T  -h+T'  r  >  1 

r 


(99b) 


Or,  putting  in  (97)  one  has 


cos  <j>  =  exp  (-i  i  I  [  (l+2r+r2 )  (E  (-1) n  r11  +  1)  ]  . 

n 


(100) 


The  square  bracketed  term,  upon  multiplying  the  factors, 
becomes 


Z 

n 


(  (_i)n  l<2±ii  rn  + 
u  }  2n+l 


(-1) 


.  ,  ,  ,  x  n+1 
n  4(n+l)r 

2n+l  ^  L) 


n  2  (n+1) 
2n+l 


n+2  . 
r  ) 


+ 


+  2r  + 


1  , 


Since  n  is  a  dummy  index  one  can  change  n  to  n-1  and  n 
to  n-2  in  the  second  and  third  terms  respectively  to 
reduce  the  sum  to  common  powers  of  r.  Thus,  one  has 
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£ 

n=l 


n  2  (n+1) 
'  2n+l 


E 

n=2 


4n  (-r)  n 
2n-l 


+  £ 
n-3 


2  (n-1 ) 
2n-3 


+  1  +  2r  + 


r 


2 


+  8  X.  (-r)n 
n=3 


1 

(2 n+1) ( 2n-l ) (2n-3) 


=  G(r) 


for  r  <  1  (101a) 


and  similarly,  for  r  >  1  one  has 


G(r) 


2,11 

3  15  r 


oo  /  i  /  \  n-1 

y  (~l/r) _ 

n^3  (2n+l)  (2n-l)  (2n-3) 


(101b) 


so  that 

cos  <J>  =  exp  [-  j  —  G(r)  I]  .  (102) 

_3 

Since  the  summation  will  be  proportional  to  n  , 
it  may  be  neglected  for  r  <1,  and  even  for  r  ~  1. 

From  a  plot  of  (de/dp)^  it  can  be  seen  (Bl  66) 
that  the  curve  may  be  approximated  by  (de/dp)  =  0.4  e 
for  the  range  1.2  <  e  <20.  Thus,  the  stopping  power 

becomes 

=  0.4  e"1^  +  k  (103) 

dp 

using  (23)  and  (73);  and  putting  these  into  (77)  and 
the  definition  for  I,  one  finds 


•- 
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o 


de 


6  =  0.4  +  ke 

e 


or 


=  -(4»i>%  b  log  (i- ^4 &?-> 


1  +  0 . 4/k e 


o 


(104) 


1  +  0.4/ke  ,  /2  ,  _ 

- -  exp  (-/-  k0) 


e_  1  +  0.4/ke 


o 


(105) 


and 


I  =  0.4  / 


de 


(0.4  +  ke )  e 


1  +  0.4/ke  -1 

lQ9  ’lToi/tr1 


(106) 


so  that 


cos  <j> 


rl  +  0.4/ke  1 

L1  +  0.4/ke  J 
o 


-krG 


(107) 


Now,  both  expressions  (105)  and  (107)  may  be  evaluated 
numerically  without  much  difficulty. 


Thus,  one  has  both  v  (or  e)  and  cos  <|>  as  ' 
functions  of  time,  0,  and  attention  can  now  be  returned 
to  the  calculation  of  F(t).  Recall  that 


CO 


f  v  ( t )  cos  <j>(t)  exp  (-t/ t )  dt 
A  & 

- =  F  (t)  =  - - 


AE 


yo 


OO 


/  v  exp(-t/x)dt 
0 


OOyK/ 

T  .  V (0) 

J  /W 

T  0  Vo 


0 


t/T 


)d0 


cos  (p  ( 0 )  exp  (- 


(108) 
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using  the  variables  0  =  t/T  and  v  =  — ^  v.  A  remark 

_  e^ 

should  be  made  here  that  v  cos  <j>  should  actually  be 


used.  However,  the  corrections  to  I  from  the  use  of 


v  cos  <J>  is  only  5%  to  10%  (B1  66)  for  r  =  1  and 

>>  (P^-)  ,  and  much  less  for  all  other  cases, 

dp  n  dp  e 


For  very  short  lifetimes  (  AE  ~  AE  and 
J  y  yo 


T<<  a  =  (2 Tic)  T/k)  a  small  quantity  F(t)  is  defined 
such  that 


F(t) 


AE 


yo 

AE 


-  AE 

X 

yo 


=  —  f  (1  -  cos  <f>)  exp  (-  — 7Ffr)d0.  (109) 

x  n  v  T/I 


Consider,  now,  the  case  where  electronic  stopping 
dominates  [(de/dp)e  >>  (de/dp)^  for  0  <  t/T].  Assume 

that 

—  exp  (— )  0.4/ke  <<  1 

r  ^  a  o 


which,  in  most  experimental  situations,  is  equivalent 
to  saying  that 


with 


0.4  = 

kso 


d  <<  1 
o 
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Thus,  from  (105)  one  finds 


1  1  -x  r  -x  .  ,  , 

e  [e  -  2d  smh  x] 

o 


d  d 


(110) 


where  x  =  QT/a,  or 


d  s;  d  e^X  [1  +  2d  eX  sinh  x] 
o  o 


(111) 


to  first  order  in  d  .  Note  also  that 

o 


d  t 

=  (-?) 
v  d 

o 


e  X  [1  -  d  eX  sinh  x] 
o 


(112) 


From  (107)  it  can  be  seen  that 


- -  r  1  +  d-, 

cos  <j)  =  Lt — ; — “ J 


-n 


1  +  d 


o 


(113) 


where  n  =  •^r—  G.  Or, 
2r 


cos  cj)~l  +  nd  -  nd  , 
Y  o 


(114) 


to  first  order  in  dQ  and  d,  so  that 


v 


-x 


-x 


v 


o 


cos  <J>  «  e  -  dQ  sinh  x  +  n  dQ  e  -  n  dQ  e 


x 


(115) 


Or, 


/V 

V 


COS  (J) 

o 


e"“9-  [1  +  |g]  sinh  (£  9) 

o 


(116) 
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Putting  this  into  expression  (108)  for  F(x)  and 
carrying  out  the  trivial  integration  one  finds 


F(t)  =  - - - ^  (l  +  jG)  - ^ - .  (117) 

1  +  x/a  £o  1  -  (x/at 


The  first  term  of  (117)  is  due  to  electronic  stopping 
and  the  second  is  due  to  nuclear  stopping,  which  itself 
is  composed  of  two  parts:  loss  of  kinetic  energy  by 
the  moving  ion  to  the  surrounding  nuclei  in  the  colli¬ 
sion,  and  scattering  (change  in  direction)  of  the 
moving  ion.  Since  G  is  of  the  order  of  unity,  one  can 
easily  see  that  the  scattering  term  dominates  for  r 
much  less  than  unity,  which  corresponds  to  the  case  of 
light  ions  slowing  down  in  a  high  mass  number  backing. 
The  above  expression  holds  for  the  range  x/a  <  0 . 3  to 
0.5. 


In  this  derivation  it  has  been  assumed  that  the 


stopping  material  is  composed  of  atoms  of  only  one  z 
and  A  ignoring  the  fact  that  many  materials  available 
are  composites  of  several  types  of  atoms.  If  771,  z^, 

A2  refer  to  the  heaviest  stopping  atom  and  Ifl^ ,  z 2i ,  A^ 

t  T- 

refer  to  the  i  lighter  atoms  then  (77)  becomes  (Bl  66) 

e 


0  =  (ify) 


de 


^  de  .  ^  ,de 


s  e2‘dt  +  l  Vaxhi 


(118) 


FI 


43 


where  0  and  e  pertain  to  the  heaviest  scattering  atom 
and 


N .  z~  .  a .  A,  + 

l  2i  l  1 _ 2 

N  z~  a  A, +  A„ . 

2  1  2i 


(119) 


de 

and  where  (^j) ^  is  the  total  specific  stopping  power 
at  an  energy  . 

For  I  one  has,  now. 


I  = 


j  A0 .  G .  j 

p  (^£.)  +  y  c  — i.  (iLf.) 

o  dp  n  .  i  A  G  'dp' ni 
r  _ - _ _ _ 


j 

e 


£t§+  f  ci  <3F>i] 


de 


(120) 


and 


d' 

n 


cos  = 


1  +  d  /d  e 
n '  e 

1  +  d  /d  e 
n'  e  o 


2r  d 


n 


(121) 


where , 


de  =  k  +  E  Ci  k±  {171^/ni) 


d  =  0.4  [1  +  £  C.  (ffl/Ws.  )  2] 
n 


<  =  °-4  t1  +  s  ci  sf  ^  {w)h 

l  2  l 


Then 


d  A„  d 1 

„  ,  x  1  n  M  2  n 

F(t)  25  1  +  x/a  d  e  (1  A,  G  d  } 


x/ql 


e  o 


n  1  -  (x/a) 


(122) 


Since  the  stopping  process  is  purely  an  atomic 
effect,  having  nothing  to  do  with  molecular  bonding, 
the  resultant  stopping  power  of  a  molecule  is  propor¬ 
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tional  to  the  sum  of  the  stopping  powers  due  to  the 


individual  atoms  of  the  molecule.  Thus,  for  a  compound 
consisting  of  molecules  of  n^  atoms  of  type  i,  the 
stopping  power  is  given  by 


(— > 

dx  molec . 


E  n.  M.  (^)  . 
^11  dx  i 

E  n.  M. 

l  l 
l 


(123) 


so,  clearly,  for  a  molecule  consisting  of  a  heavy  and 

light  atom,  the  greatest  contribution  to  the  stopping 

will  come  from  the  heavy  atom. 

It  has  been  shown  previously  that  where  large 

dn 

angle  scattering  is  negligible  (-5—  <<  1,  x/a  <  1) 

aee 

electronic  stopping  dominates  and  that  this  should  be 
proportional  to  the  velocity  of  the  ion*. 


/dE >  _  dE  _  , 

Kdx}e~  dx 


v  (t) 


(124) 


where  k  is  a  constant  of  proportionality 


dE 

dx 


One  notes  that 


=  (ib/d£)  =  M1  v<^>/v  =  M 


dv, 


dv  (t) 


ldt  dt 


dt' 


1  dt 


(125) 


Therefore  by  combining  (125)  and  (124) ,  and  integrating 


one  has  that 


•- 


t) 


(126) 


v  =  v  exp  (-  — 
o  ^  M1 

where  vq  is  the  initial  speed  of  the  ion. 
Define  a  parameter  a  such  that 
M-^/a  ^  k 

or,  in  dimensionless  variables 


a  =  (2  M)^  T/k 


which  is  known  as  the  characteristic  electronic  slowing 
down  time . 

Ignoring  the  scattering  term  cos  <j>’  in  (17)  one 


has 


F(t)  =  /  v  exp  (-t/a)  exp(-t/x)  dt 

o  0 

1 

1  +  T/a 


(127) 


for  velocity  proportional  stoppers. 

Robertson  (Ro  69)  has  shown  that  for  the  parti- 
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cular  case  of  K  stopping  m  K  I,  the  correction  to 
F  (t)  due  to  strong  nuclear  stopping  is  about  25%. 

Since  the  masses  of  both  the  recoil  and  stopping  nuclei 
are  the  same  one  would  expect  nuclear  stopping  to  be 
large  for  this  case  and,  hence,  the  25%  error  in  (127) 
for  other  cases  to  be  an  upper  limit. 


■ 
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CHAPTER  2 

DETECTORS 

2 . 1  Neutron  Detector 

The  neutron  detector  consists  of  a  liquid  organic 
scintillator  (NE  218)  in  a  chamber  viewed  by  an  RCA  8575 
photomultiplier  providing  a  fast  50  ohm  impedance  dynode 
output  and  a  slow  linear  100  ohm  anode  output.  Inputs 
to  this  unit  are  the  high  voltage  connection  to  the 
photomultiplier  section  and  the  pre-amplifier  and  dis¬ 
criminator  inputs.  The  fast  signal,  so  called  because 
it  is  taken  from  a  dynode  in  the  photomultiplier  ampli¬ 
fication  chain  and,  thus,  is  spared  R-C  shaping  involving 
long  time  constants,  is  used  for  timing  purposes  after 
it  is  passed  through  a  walk-free  timing  discriminator 
(rise  time  of  this  timing  signal  T  <  2  n.s.  and  the 
pulse  width  T^  ~  10  n.s.).  The  anode  signal  is  taken 
to  a  preamplifier  and  integrated  so  that  the  output  rise 
time  is  approximately  50  n.s.  and  is  suitable  for  pro¬ 
cessing  by  an  R-C  shaping  amplifier. 

The  scintillator  responds  to  both  gamma  rays 
and  neutrons  so  that  for  n-y  coincidence  work,  methods 
must  be  found  to  identify  the  gamma  rays  and  also,  if 
possible,  to  separate  the  gamma  ray  linear  signals  from 
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the  neutron  signals  (pulse  shape  discrimination) .  The 
pulse  shape  discrimination  determines  what  radiation 
dependent  effects  combine  to  produce  the  scintillations 
in  order  to  correlate  the  observed  events  with  the  type 
of  radiation  incident  on  the  detector. 

Organic  scintillators  are  composed  of  hydrogen 
(proton)  rich  molecules  (the  specific  case  here,  NE  218, 
has  a  proton-carbon  ratio  of  1.4).  Scintillations  are 
produced  when  either  a  proton  or  an  electron  produces 
"excited"  molecules  along  its  flight  path  through  the 
liquid.  This  "excitation"  is  of  two  types  (Br  59)  - 
electronic  excitation  and  ionization. 

Electronically  excited  molecules,  M* ,  decay  by 
two  processes: 

1.  M*  -*  M  +  photon 

2.  M*  ->  M  +  heat  (vibrational  energy)  (A) 

Ionization,  however,  involves  an  additional  .step: 

M+  +  e""  -*  M*  (B) 

with  M*  decaying  by  either  process  1  or  2  of  (A) . 

Process  (B)  involves  the  recombination  of  the  ion 
with  an  electron  and  thus  is  governed  by  the  relatively 

-7 

slow  diffusion  of  electrons  (recombination  time  >  10 
secs.) ,  whereas,  process  (A)  takes  of  the  order  of 
10"9  secs,  to  complete.  The  light  intensity  of  both 
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processes (A)  and  (B)  are  exponential  decays  in  time. 

It  has  been  shown  (Ta  51)  that  the  light  output 
or  pulse  height  in  an  organic  scintillator  is  not 
proportional  to  the  particle  energy  and,  in  fact,  is 
less  for  a  proton  than  for  an  electron  of  the  same 
energy.  The  light  output  is  actually  a  function  of 
the  density  of  M  and  M  along  the  trajectory  of  the 
particle  and  it  is  here  that  different  modes  of  quench¬ 
ing  (other  than  1  of  (A) ) ,  through  interactions  (of  type 
2,  for  example)  between  excited  ions  of  molecules,  can 
occur,  thus  lowering  the  light  output  since  type  1  is 
the  only  photon  producer.  Clearly,  a  heavier  particle 
is  going  to  leave  a  higher  density  of  excited  molecules 
in  its  wake  than  a  lighter  particle  of  the  same  energy 
and  the  light  output  from  protons  is  less  than  that  from 
electrons  since  the  former  case  the  greater  density  of 
excited  molecules  permits  quenching  by  processes  other 
than  type  1. 

The  scintillations  may  now  be  divided  into  two 
types:  a  fast  component,  F,  which  is  attributable  to 
process  1  of  (A),  and  a  slow  component,  S,  from  process 
(B) .  The  decay  time  of  F  is  usually  several  nanoseconds 
while  for  S  it  is  about  100  n.s.  Further,  as  has  been 
stated  before,  one  expects  more  photons,  f  ,  from  the 
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electrons  than  the  number,  f  ,  from  the  protons  (f  <  f  ) . 

p  p  e 

The  production  of  the  slow  component,  S,  has  two 
effects:  one  is  that  it  is,  by  the  nature  of  the  recom¬ 
bination  process,  slow,  and  the  other  is  because  it  is 
so  slow,  the  M*  formed  through  process  (B)  arrive  in 
a  situation  where  the  density  of  initially  excited 
molecules  (usually  formed  through  (A)  and  hence  have 
fast  decay  times)  is  much  lower  than  it  was  at  the 
beginning  of  excitation  so  that  the  probability  of  type 
1  occurring  is  greatly  enhanced  over  process  2.  There¬ 
fore,  by  the  previous  argument,  more  "slow"  photons  S’ 
are  usually  expected  from  protons  than  the  number  from 
electrons.  Although  it  is  still  possible  for  <  Se 

the  ratio  S/S  >f/f  so  that  S  / (f  +  S  ) ,  the  pro- 
p  e  p  e  P  P  P 

portion  of  S  from  proton  scintillations  is  greater  than 

S  / ( f  +  S  )  from  electron  scintillations. 
e/  e  e 

Now,  one  need  only  recall  that  gamma  rays  produce 
recoil  electrons  (photoelectrons,  Compton  electrons,  or 
pair-produced  electrons  and  positrons)  and  neutrons 
produce  recoil  protons,  to  conclude  that  gamma-originated 
scintillations  will  have  less  slow  component  than  neutron- 
originated  scintillations  (figure  4) . 

In  delay  line  shaping  the  input  signal  is  divided 
into  two  branches:  in  one  branch  the  signal  current  is 
unaffected  but  in  the  other  the  inserted  delay  line 
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Figure  4.  All  vertical  scales  are  pulse 
amplitudes  and  all  horizontal  scales  are 
time . 

(a)  Solid  line  represents  a  neutron  scin¬ 
tillation  with  an  exponential  decay 
having  fast  (F)  and  slow  (s)  com¬ 
ponents  and  the  dotted  continuation 
completes  the  gamma  ray  produced 
scintillation  (almost  entirely  fast 
component) . 

(b)  Preamplification  of  (a) . 

(c)  ,  (d) ,  and  (e) ,  (f )  show  the  process 

of  single  delay  line  shaping  applied 
twice  resulting  in  separate  crossover 
points  (tp,  t  )  for  the  gamma  ray 
originated  pulse  and  the  neutron  pulse. 


<  2  O _ i 


51 


inverts  and  delays  the  signal.  Resistors  are  placed 
in  both  branches  to  compensate  for  the  attenuation  of 
the  delay  line.  The  currents  in  the  branches  are  then 
summed  at  a  common  junction,  producing  a  shaped  output 
signal.  Double  delay  line  shaping,  of  course,  performs 
the  process  twice.  The  effect  is  shown  (figure  4)  for 
two  input  pulses  (one  with  a  large  slow  component  and 
one  with  very  little  slow  component) .  The  effect  of  the 
slow  component  is  to  cause  an  incline  to  the  plateau 
region  of  the  signal.  After  double  delay  line  shaping 
this  results  in  a  displacement,  in  time,  of  the  zero 
crossover  of  the  signal  and  by  using  a  zero  crossing 
discriminator,  fast  timing  signals  will  be  generated  at 
the  two  crossover  points.  Neutron  identification  then 
proceeds  in  the  following  manner. 

A  zero  time  mark  may  be  derived  from  the  neutron 
detector  discriminator  output  which  takes  its  signal 
from  a  dynode  at  the  onset  of  charge  collection  at  the 
anode  (it  is  the  same  for  either  the  neutron  or  the 
gamma  ray)  and  is  used  to  start  a  time-to-pulse  height 
converter  (TAC) .  The  timing  signals  derived  from  the 
crossover  of  the  shaped  linear  pulse  stop  the  TAC  and 
the  time  spectrum  obtained  will  contain  two  peaks  -  one 
for  neutrons  and  one  for  gamma  rays,  and  may  be  used  as 
a  neutron  identifier  for  the  detector. 
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The  scintillator-photomultiplier  is  used  most 
advantageously  in  fast  timing  applications  because  of 
its  excellent  timing  resolution.  This  resolution, 
however,  is  governed  by  four  main  factors  (Re  61) : 

a)  variation  of  the  time  of  interaction  of  the 
radiation  with  the  scintillator  and  the 
amount  of  energy  deposited  therein  (detected 
thickness  contributes  to  this) ; 

b)  the  finite  decay  times  of  the  light  emitting 
states  in  the  scintillator  and  variation  of 
the  times  of  photon  arrival  at  the  photo¬ 
multiplier  cathode  due  to  optical  reflections 
and  absorption  within  the  scintillator  chamber 
which  are  determined  by  the  geometry  of  the 
chamber  and  by  the  optical  coupling  of  the 
scintillator  to  the  photomultiplier; 

c)  the  variation  of  the  transit  time  of  the  photo¬ 
electrons  in  the  photomultiplier  due  to  the 
different  path  lengths  and  the  variation  of 
the  critical  energy  and  angle  of  the  secondary 
electrons;  and 

d)  the  jitter  and  uncertainty  of  the  times  of 
triggering  of  associated  electronics. 
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2.2  Lithium  Drifted  Germanium  Detector 


A  simple  semiconductor  detector  consists  of 
a  junction  between  n-type  material,  where  virtually 
all  the  charge  carriers  are  electrons,  and  p-type 
material,  where  they  are  all  holes  (this  is  basically 
a  diode) .  Where  the  two  are  joined  some  of  the  holes 
from  the  p-layer  diffuse  into  the  n-layer,  and  some 
of  the  electrons  from  the  n-layer  diffuse  into  the 
p-layer,  the  process  constituting  an  electric  current 
flowing  from  the  p  to  the  n  region  and  clearly  cannot 
continue  indefinitely  without  power.  The  current  is 
inhibited  by  the  built-in  voltage  arising  from  ionized 
impurities  in  each  region.  If  this  diode  is  reverse- 
biased  the  holes  will  be  forced  into  the  p-region  and 
electrons  into  the  n-region  creating  a  depletion  layer 
which  is  charge  free.  One  should  note  that  one  way  of 
creating  intrinsic  material  from  substances  such  as 
germanium  or  silicon  having  p-type  impurities  is  to 
drift  the  material  with  an  n-type  material  such  as 
lithium. 

To  find  the  depletion  layer  thickness  one  has 
to  solve  Poisson's  equation.  For  p-type  material  this 
becomes  (in  one  dimension)  (De  63) 
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where  V (x)  is  the  potential  difference  across  the  p- 

region,  x  is  the  dielectric  constant  of  the  crystal, 

N  is  the  acceptor  impurity  concentration,  n  and  p  are 

the  local  electron  and  hole  densities. 

Since  in  actual  fact  n  and  p  are  not  constants, 

solution  of  (1)  will  present  a  rather  tedious  problem, 

so  it  will  be  assumed  that  n  and  p  are  small  compared 

to  N  throughout  the  depletion  region  (a  poor  appro- 
a 

ximation) .  Integrating  (1)  twice  one  finds  that 

4  t r  e  N 

V  (x)  =  - - — -  (x2-2x  xp)  +  Vq  (2) 


where  x  is  the  termination  of  the  depletion  region 
P 

where  N  +  n  -  p  =  0  and  VQ  =  V(0)  is  the  potential 
at  x  =  0,  the  junction  interface.  Thus  at  x  =  x  , 

V  -  VQ  =  V  ,  the  potential  difference  across  the  p- 
region.  Therefore, 


V  X 

P  * _ 

2tt  e 
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is  the  thickness  of  the  p-region  and. 
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is  the  thickness  of  the  n-region,  where  NQ  is  the  donor 
impurity  concentration  and  is  the  potential  across 
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the  n-layer.  In  actual  cases,  the  ionized  impurity 
concentration  is  very  different  on  the  two  sides  of 
the  junction  and  the  depletion  thickness  resides 
almost  entirely  within  one  type  of  material.  For 
n-type  material  the  width  will  be  given  by 


_ X__  X 

2  it  e  Nq 


(5) 


The  resistivity  of  the  n-type  material  will  be 


given  by 


P  = 


N„  e  p 
D  ^n 


(6) 


where  u  is  the  electron  mobility  in  the  crystal  and, 
Kn 


U 


n 


(7) 


E  being  the  applied  electric  field  and  vg  the  velocity 
of  the  electron,  and  is  given  by 


E  =  J  p  ,  (8) 

J,  of  course,  is  the  current  density  and  may  be  expressed 
in  terms  of  NQ  and  ve: 


J 


(9) 
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Equations  (7),  (8),  and  (9)  combine  to  give  (6),  and 

(6)  with  (5)  gives 


x 


(10) 


This  brief  description,  of  course,  is  for  a  planar 
detector ,  whereas  the  actual  detector  used  was  co¬ 
axial,  the  situation  being  considerably  more  complicated. 
The  principle,  however,  is  basically  the  same  and  one 
usually  has  a  p-type  core  with  lithium  drifted  radially 
inwards  to  form  an  n-type  cylindrical  layer. 


The  depletion  region  behaves  in  an  analogous 


manner  to  an  ionization  chamber  since  a  charged  particle 
incident  upon  the  region  loses  energy  creating  hole- 
electron  pairs  which  are  swept  out  by  the  applied  field, 
the  charge  being  collected  as  a  function  of  time.  On 
the  average,  it  takes  about  3.5  eV  of  energy  lost  by 
the  incident  particle  to  create  one  hole-electron  -pair 
in  silicon  and  since  this  is  constant,  the  energy  of  the 
particle  will  be  proportional  to  the  total  charge  it 
produces.  Photons  are  detected  in  the  crystal  by  the 
transfer  of  their  energy  to  charged  particles  in  the 
detector.  This  interaction  is  governed  by  three 
processes  by  which  photons  interact  with  matter:  the 
photoelectric  effect,  where  all  the  energy  of  the  photon 
is  transferred  to  a  photoelectron;  the  Compton  effect. 


Figure  5.  The  basic  time  walk  problem 


of  a  simple  discriminator  is  illustrated 
for  two  pulses  of  amplitudes  and  V ^ • 
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where  part  of  the  photon  energy  is  transferred  to  an 
electron  and  the  photon  of  reduced  energy  is  scattered; 
pair  production,  where  the  photon  energy  is  used  to 
create  an  electron-positron  pair  having  some  kinetic 
energy.  Thus,  the  Ge  (Li)  detector  may  be  used  as  an 
energy  sensitive  device  for  the  detection  of  photons 
and  charged  particles. 

It  is  possible  to  use  pulses  from  Ge  (Li) 
detectors  for  timing  purposes,  as  well,  providing  the 
signals  are  processed  in  the  proper  manner  to  eliminate 
time  walk  problems.  This  problem  is  typical  of  leading- 
edge  timing  systems  where  a  timing  signal  is  produced 
when  the  input  signal  amplitude  rises  above  a  pre-set 
discriminator  level.  Variations  in  amplitude  lead  to 
shifting  of  the  timing  point  along  the  time  axis 
(figure  5) .  The  problem  can  be  very  drastically  reduced 
if  one  could  cause  the  discriminator  level  to  move  up 
and  down  with  the  pulse  amplitude  so  that  the  time  mark 
is  generated  at  a  constant  fraction  of  pulse  height 
(constant  fraction  of  pulse  height  trigger  or  C.F.P.H.T.). 

Another  problem  that  spoils  the  time  resolution 
of  the  Ge  (Li)  detector  is  the  variation  of  pulse  shapes 
due  to  the  fact  that  charged  particle  interactions  occur 
at  different  locations  within  the  depletion  region 


Figure  6 .  A  planar  detector  of  width 

W  separating  the  n+  and  p  layers 

(depletion  region)  is  shown  biased 

with  V  volts.  x  is  the  location  of 

the  ionization  and  the  motion  of  the 

+ 

resulting  holes  (h  )  and  electron 
(e  )  is  indicated. 

The  pulse  shapes  for  ionizations 
occurring  at  positions  A  and  B  in  the 
depletion  layer  are  depicted  for  the 
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(Sa  68).  To  see  this  more  clearly  (figure  6a)  consider 
a  planar  detector  with  a  depletion  region,  i,  of  total 
width,  W,  and  with  a  charged  particle  interaction 
occurring  at  a  distance  xq  from  the  n  layer  producing 
nQ  hole-electron  pairs  each  of  charge  q  and  mobility 
y^  and  yg  respectively.  t  is  the  time  after  ionization 
and  the  quantities  tQ  and  t^  will  be  defined  as 


follows , 


t  =  t 


W  x 


for  t  ^ 
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he  v 


t  = 


W  x 
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y  v 
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t,  =  t 
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^h  V 
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for  t  > 


for  t< 
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since  the  velocity  for  a  charge  carrier  of  mobility  y 
in  an  electric  field  E  =  V/W  is  given  by 


v  =  y  E 


U  V 

W 


(ID 


so  that  the  time,  t,  to  travel  a  distance  x  is 

I 

=  x  =  wx  _  (12) 

v  y  V 


Let  T  be  the  time  to  travel  over  the  entire  width 


W.  Therefore, 
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(13) 


for  electrons  and. 


W‘ 


^  v 


(14) 


for  holes.  Hence  the  ratio  of  the  charge  collected 
after  time,  t,  to  the  total  charge  produced  at  ioniza¬ 
tion  is  given  by 


Q  (t) 

n  q 
o  ^ 


(15) 


a  plot  of  which  is  seen  in  figure  6b  where  some  typical 
pulse  shapes  are  illustrated.  The  pulse  shapes  are 
characterised  by  an  initial  fast  rise  until  one  type 
of  carrier  has  been  swept  out  of  the  depletion  region, 
followed  by  a  rise  at  half  of  the  former  rate  of  charge 
collection  until  the  other  type  of  carrier  reaches  its 
appropriate  layer.  Of  course,  the  possibility  of' 

Compton  scattering  and  multiple  interactions  has  not 
been  considered,  the  effect  being  to  complicate  the 
shapes  of  the  pulses.  There  will  still  be  a  fast  and 
a  slow  component  to  the  pulses,  however. 

For  the  case  of  a  co-axial  detector  having  an 
n+  layer  and  a  p-type  core  the  situation  is  made  quite 
complicated  because  of  the  1/r  dependence  of  the  electric 
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field.  Electrons  are  accelerated  outwards  by  a 
decreasing  field  and  holes  inwards  by  an  increasing 
field.  The  result  is  that  ionizations  occurring 
close  to  the  n+  layer  tend  to  produce  pulses  with 
slow  rise  times  whereas  events  midway  between  the  n+ 
and  p  layers,  in  a  considerably  wide  band,  produce 
pulses  with  relatively  fast  rise  times.  Events  close 
to  the  core  produce  pulses  with  initially  fast  rise 
times  to  about  Q(t)/nQq  z  0.5  with  slow  components 
taking  over  until  the  total  charge  is  collected. 

The  effect  of  this  pulse  shape  distribution  is 
to  cause  non-amplitude  related  timing  uncertainties 
which  cannot  be  entirely  compensated  for  by  the 
C.F.P.H.T.  Placing  an  annulus  of  lead  over  the  outer 
regions  of  the  detector  will  help  suppress  slow  rising 
pulses  by  partial  elimination  of  ionizations  in  that 
location,  but  the  effect  will  not  be  great  and  the 
disadvantage  will  be  a  drastic  reduction  in  the  Ge  (Li) 
efficiency  since  from  purely  geometrical  considerations, 
one  can  see  that  a  large  fraction  of  the  total  number 
of  events  originate  in  this  region.  Another  method  to 
improve  the  time  resolution  might  be  to  use  a  rise  time 
spectrometer  and  a  computer  to  select  pulses  according 
to  their  rise  times  so  that  separate  time  pickoffs  can 
be  used  for  groups  of  pulses  with  similar  rise  times. 
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CHAPTER  3 

EXPERIMENTAL  ASPECTS  OF  THE  53Mn  LIFETIME 

MEASUREMENTS 

3.1  Introduction 

This  experiment  is  concerned  with  the  measurement 

5  3 

of  the  lifetimes  of  the  bound  levels  of  Mn  through  the 
5?  53 

Cr  (p,  n  y)  Mn  reaction  using  the  DSAM  method. 
Although  the  DSAM  method  may  be  employed  without  the 
neutron-gamma  ray  coincidence,  the  inclusion  of  this 
condition  greatly  simplifies  the  analysis  of  the  data. 
This  results  from  the  reaction  kinematics  which  show 
that  if  only  a  narrow  cone  of  neutrons  is  detected  in, 
say,  the  backward  direction  (from  the  incident  beam 
direction)  then  the  de-excitation  gamma  rays  in  coin¬ 
cidence  with  these  neutrons  will  be  from  recoil  nuclei 
travelling  in  a  much  narrower  cone  in  the  forward 
direction  (for  this  experiment  with  5.5  MeV  protons 
the  recoil  nuclei  will  subtend  a  half-angle  of  less 

than  4°  if  the  neutron  detector  subtends  a  half-angle 

■% 

of  20°) .  Thus,  the  effect  on  the  gamma  ray  Doppler 
shifts  due  to  the  initial  angular  dispersion  of  the 
recoil  nuclei  is  much  more  easily  calculated  for  the 
case  of  the  narrow  cone  than  for  an  isotropic  angular 


Figure  7.  The  electronics  set  up  is  shown. 

(1)  R-C  amplifier  (input  from  the  Ge  (Li) ) 
(Ortec  Model  450) . 

(2)  Single  channel  analyzer 
(Ortec  Model  406) . 

(3)  Timing  filter  amplifier  (accepts 
timing  signals  from  the  Ge  (Li) ) 

(Ortec  Model  454) . 

(4)  Constant  fraction  timing  discriminator 
(C.F.P.H.T.)  (Ortec  Model  453). 

(5)  Gate  and  delay  generator 
(Ortec  Model  416). 

(6)  Time  to  amplitude  converter 
(Ortec  Model  437) . 

(7)  Double  delay  line  amplifier 
(Ortec  Model  410) . 

(8)  Fast  zero  crossing  discriminator. 

(9)  Time  pickoff  control 
(Ortec  Model  403A) . 

(10)  Universal  coincidence  units 
(Ortec  Model  418) . 

(11)  Timing  single  channel  analyzer 
(Canberra  Model  1436) . 
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distribution  (no  neutron-gamma  ray  coincidence) . 

Another  advantage  of  the  coincidence  method  is  that 
since  the  only  recoil  nuclei  considered  are  those 
confined  to  the  forward  direction  they  will  have  the 
highest  initial  velocity  and  will,  therefore,  result 
in  the  greatest  possible  Doppler  shifts. 

In  this  chapter  the  electronics  system  required 
to  obtain  a  calibrated,  neutron-coincident  gamma  ray 
spectrum  (using  pulse  shape  discrimination  in  the 
neutron  detector  as  discussed  in  Chapter  2)  will  be 
described  along  with  the  mechanical  arrangement  of 
the  detectors  and  the  type  of  targets  used  in  the 
experiment . 

3.2  Electronics 

The  electronics  system  is  shown  in  figure  7. 

The  gamma  ray  pulses  from  the  Ge  (Li)  detector  are 
put  into  an  R-C  amplifier  with  a  shaping  constant  of 
2  y . s . ,  the  outputs  going  to  a  single  channel  analyzer 
(SCA)  and  to  the  analog-to-digital  converter  (ADC) . 

A’  window  is  set  on  the  desired  gamma  spectrum  with 
the  SCA  the  output  of  which  goes  to  a  coincidence  unit. 
The  gamma  pulses  from  the  detector  are  also  put  into  a 
timing  and  filter  amplifier  (TFA)  to  be  properly  R-C 
shaped  for  input  to  the  CFPHT .  The  slow  logic  signal 
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thus  derived,  is  then  passed  on  to  two  gate  and  delay 
units,  the  delayed  marker  of  one  going  to  the  gamma- 
gamma  TAC ,  the  other  going  to  the  start  input  of  the 
n-y  TAC  the  stop  coming  from  the  fast  neutron  detector 
output.  The  output  of  this  TAC  gives  the  n-y  time 
spectrum.  The  extra  gate  and  delay  in  the  fast  neutron 
side  is  to  provide  maximum  flexibility  in  the  adjust¬ 
ment  of  the  delays.  The  gamma  ray  pulses  from  the 
sodium  iodide  (Na  I)  detector  (which  is  shielded  with 

lead  so  that  it  views  only  a  ^Co  source  whereas  the 

53 

Ge  (Li)  detector  views  both  the  Cr  target  and  the 
source)  are  double-delay-line  shaped  and  then  a  timing 
SCA  (TSCA)  is  used  to  put  a  window  on  the  two  source 
lines.  The  fast  output  of  the  TSCA  goes  to  stop  the 
y-y  TAC,  the  output  of  which  goes  to  the  same  coincidence 
unit  as  before  (coincidence  requirement  of  2).  Hence, 
a  gamma  ray  coming  from  the  source  (within  the  6^Co 
window)  and  a  gamma  from  the  target  (within  the  spectrum 
window)  will  cause  the  coincidence  unit  to  deliver  an 
output  going  to  the  final  coincidence  unit  and  to  the 
2"  route  of  the  kicksorter,  meaning  that  the  gamma  spectrum 
with  the  source  spectrum  included  will  be  put  into  the 
next  2048  channels  from  the  primary  n-y  coincidence 
spectrum.  The  neutron  linear  signal  is  passed  through 
a  fast  zero  crossing  discriminator  (FZCD)  which  delivers 
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early  pulses  for  gamma  ray  pulses  and  late  pulses  for 
neutron  signals  into  the  stop  input  of  the  neutron 
identifier  TAC  which  is  started  by  the  neutron  detector 
discriminator  pulse  occurring  for  both  neutrons  and 
gamma  rays.  The  output  of  the  TAC  will  contain  two 
peaks  -  one  for  neutrons  and  the  other  for  gamma  rays, 
counted  in  the  detector.  Windows  are  placed  on  the 
neutron  identifier  peak  and  on  the  neutron  peaks  of 
the  n-y  time  spectrum  using  SCA's,  the  pulses  of  which 
are  placed  in  coincidence  with  the  gamma  spectrum 
window  (3-fold  coincidence  requirement)  in  a  coincidence 
unit  which  provides  a  true  n-y signal.  This  signal  and 
the  y-y  signal  are  both  fed  into  a  common  gate  having 
a  single  coincidence  requirement  so  that  either  signal 
opens  the  ADC  gate.  The  gamma  rays  which  are  in  coin¬ 
cidence  with  the  neutrons  appear  in  the  first  2048 
channels  of  the  kicksorter.  A  feature  of  the  n-y  TAC 
that  provides  a  true  start  has  been  used.  The  true 
start  signal  is  delayed  and  then  used  to  strobe  all 
three  TAC ' s  at  the  same  time,  if  and  only  if  a  stop 
signal  at  the  n-y  TAC  has  been  preceded  by  a  legitimate 
start  signal.  This  feature  is  advantageous  where  high 
count  rates  are  prevalent  such  as  the  neutron  detector 
since  it  detects  both  neutrons  and  gamma  rays.  Record¬ 
ing  the  independent  y-y  spectrum  with  its  unchanging 
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^Co  peaks  allows  one  to  monitor  any  gain  shifts  that 
may  occur  as  a  result  of  electronic  instabilities  in 
the  system. 

3.3  Experimental  Arrangement 

The  neutron  detector  was  placed  at  approximately 
180°  and  15  cm.  back  from  the  target  (a  hole  in  the 
center  of  the  detector  permitted  this)  and  the  entire 
face  of  the  detector  opposite  the  target  was  shielded 
with  lead  to  reduce  the  gamma  ray  flux  into  the 
scintillator.  The  fast  and  linear  outputs  were 
connected,  with  50  ft  and  100  ft  coaxial  cable,  respec¬ 
tively,  to  the  electronics  in  the  control  room.  The 
detector  was  then  biased  at  -2250  volts  and  the  dis¬ 
criminator  was  set  just  above  the  self-tr iggering  level. 
The  Ge  (Li)  detector  was  placed  on  a  rotatable  platform, 
the  center  of  rotation  lying  directly  below  the  target. 
Lead  was  placed  over  the  face  of  the  detector  to  reduce 
the  excessively  high  count  rate  produced  by  low  energy 
gamma  rays,  and  the  detector  was  biased  after  all  the 
preamplifier  connections  had  been  made. 

About  30  cm.  directly  above  the  Ge  (Li)  detector 
was  placed  the  Na  I  detector,  the  ^Co  source  being 
situated  midway  between  the  two.  The  Na  I  detector 
was  shielded  with  lead  to  prevent  viewing  of  the  target- 
originated  gamma  rays.  The  system  permits  spectra  to  be 
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taken  at  0°  and  ~  120°  to  obtain  shifts  of  the  order 

of  1.5Ae  (see  figure  8  and  figure  9). 

3 . 4  Targets 

Two  types  of  targets  were  made.  One  type  involved 

.  2 

a  simple  evaporation  of  a  small  amount  (5  0.5  mg/cm  ) 

53 

of  Cr  enriched  to  96%  Cr  onto  a  support  backing  of 

tantalum  (0.005"  thick).  The  range  of  the  recoil  ions 
53 

in  Cr  (see  Chapter  4)  was  determined  to  be  approxi- 

2 

mately  0.01  mg/cm  for  an  initial  recoil  energy  of  about 

300  KeV.  Also  it  was  calculated  that  the  beam  would 

lose  about  5  KeV  of  its  original  energy  in  passing 

through  the  target,  giving  very  little  uncertainty  in 

the  recoil  ion  starting  velocity.  However,  in  order 

to  be  absolutely  sure  all  the  ions  are  stopped  in  the 

target  material  a  target  thickness  of  several  hundred 

times  the  recoil  range  is  desirable. 

Thus  a  second  type  of  target  was  made  by  first 

evaporating  a  layer  of  natural  chromium  to  a  thickness 

2 

of  approximately  6  mg/cm  (natural  chromium  contains 

53  2 

about  6.5%  Cr)  followed  by  about  0.5  mg/cm  of 

5  3 

Cr.  The  change  in  slowing  down  time  caused  by  the 

differing  mass  numbers  will  effectively  be  far  less 

than  2%  since  the  great  majority  of  the  ions  will  be 

53 

stopped  in  the  Cr  layer.  The  beam  dispersion,  however. 
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will  be  approximately  50  KeV  resulting  in  a  small  un¬ 
certainty  in  the  starting  velocity,  and  may  necessitate 
interpolation  between  F(t)  curves  resulting  from  diffe¬ 
rent  assumed  starting  velocities. 


Figure  8.  The  side  of  the  basic  mechanical 
arrangement  of  the  experiment  is  shown: 

(1)  portion  of  the  Na  I  detector  shielded 
with  lead  from  the  target  gamma  rays; 

ft  0 

(2)  uCo  source  location; 

(3)  portion  of  the  Ge  (Li)  detector; 

(4)  beam  line; 

(5)  neutron  detector  shielded  on  both  faces 
with  lead  (hatched  line) ; 

(6)  target  inside  the  target  pot  assembly 
(joined  to  the  beam  line  behind  the 
neutron  detector) . 

The  entire  assembly  of  (1) ,  (2) ,  and  (3) 

rotates  together  about  the  target  (dotted 
line  indicates  the  vertical  axis  of  rotation) . 
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Figure  9.  The  top  view  of  the  experi¬ 
mental  assembly  is  shown. 


(1)  Ge(Li)  and  Na  I  detectors; 

(2)  target  and  target  pot; 

(3)  neutron  detector  with  photo¬ 
multiplier  , 

(4)  attached; 

(5)  beam  line. 

An  outline  of  the  correlation 
table  is  illustrated. 
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CHAPTER  4 

DATA  ANALYSIS 

4.1  Range  and  Stopping  Powers 

A  program  (DEDX)  was  written  in  APL  to  determine 
the  stopping  power  and  range  according  to  the  Lindhard 
or  Firsov  formulae.  The  range  calculations  for  the 
recoil  ions  in  the  target  material  are  necessary  in 
order  to  determine  the  target  thickness  required  to 
stop  all  the  ions  in  the  substance.  A  target  that 
is  too  thin  will  lead  to  some  vacuum  recoils  which 
make  the  DSAM  results  incorrect.  On  the  other  hand, 
a  target  that  is  too  thick  will  cause  excessive  beam 
energy  dispersion  which  leads  to  uncertainty  in 
determining  the  starting  velocity  of  the  recoils. 

The  accuracy  of  the  information  from  these  formulae 
is  better  than  20%. 

4.2  Reaction  Kinematics  and  F(t) 

A  Fortran  program  (DSAM  1A)  calculating  rela¬ 
tivistic  kinematics  was  used  to  determine  the  initial 
velocity  of  the  recoil  ions.  The  neutron  detector 
subtends  a  half-angle  of  about  20°  so  that  only  neutrons 
falling  within  this  were  considered  which,  of  course. 
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fixed  the  direction  and  energy  of  the  recoil  ions  for 

a  particular  excited  gamma  ray  level.  Since  the 

neutron  angles  considered  in  the  kinematics  were  180°, 

170°,  and  160°,  an  average  of  the  three  corresponding 

recoil  starting  velocities  was  taken.  The  variation 

among  these  velocities  for  a  particular  level  were 

typically  less  than  3%,  this  fact  being  the  empirical 

justification  for  the  simple  averaging  process.  These 

were  weighted  with  the  cos  0  term  to  determine  the 

full  kinematic  shift  for  0  equal  to  0°  and  120°,  and 

0°  and  90°  (AE  ) . 

yo 

E^q,  of  course,  is  the  unshifted  gamma  ray 
energy  and  can  be  determined  accurately  by  measure¬ 
ment  of  the  gamma  ray  at  a  forward  and  a  backward  angle 
At  the  forward  angle,  0^,  the  energy,  E  will  be 
Doppler  shifted  up  from  E^q.  Thus  one  has 


E  £  =  E  (1  +  A  cos  0.)  CD 

yf  yo  f 

where  A  is  the  Doppler  shift  attenuation  factor. 
Similarly  for  the  backward  measurement  one  has 


Hence 


E 


yb 


E  (1  -  A  cos  0,  ) 
yo  b 


(2) 


E_e  +  E,=2E  +E  A  (cos  Q 
yf  yb  yo  yo  f 


cos  0^) 


(3) 
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and 


Eyf  "  Eyb  =  EyO  A  (COS  9f+  C°S  °b)  '  (4) 

Therefore,  using  (3)  and  (4),  the  unshifted  gamma  ray 
energy  will  be  given  by  (Ro  69) , 


E 

Y° 


E  _  +  E  , 
yf  yb 


cos 

^cos 


0^-  COS  0 
0£+  cos  0 


~)  .  (5) 

b 


The  F(t)  values  were  calculated  using  an  ALGOL 
program  written  by  Robertson  (Ro  69)  using  various 
approximations  (see  Chapter  1,  nuclear  stopping)  and 
truncations  of  integrals.  These  values  agree  with 
Blaugrunds  calculations  to  within  10%  over  most  of 
the  range  of  F(x).  The  values  were  plotted  on  semi- 
logarithmic  graph  paper  as  a  function  of  t  (figure  15) 
and  are  used  to  determine  the  lifetimes  of  the  states 
studied.  The  accuracy  of  the  calculated  F(t)  values 
should  be  within  15%  to  20%. 

4.3  Centroid  Analysis 

An  APL  program  (BACK)  was  written  to  manipulate 
single  gamma  ray  peaks  (background  included)  and  per¬ 
form  background  subtraction  (two  independent  linear 
background  fits  may  be  chosen,  one  for  the  high  energy 
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side  of  the  peak  and  one  for  the  low  energy  side,  these 
being  used  to  help  determine  manually  chosen  points 
through  which  a  linear  fit  is  made  to  form  the  background 
under  the  peak) ,  peak  plotting  (with  errors) ,  peak  summa¬ 
tion  with  errors,  and  centroid  calculation  with  its 
inherent  error.  The  contribution  to  the  centroid  error 
by  uncertainties  in  determining  the  proper  background 
was  estimated  by  variation  of  the  background  within 
reasonable  limits  and  examination  of  the  resultant 
change  in  the  centroid.  The  total  error,  then,  is  simply 
the  sum  in  quadrature  of  these  errors. 
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CHAPTER  5 

THE  53Mn  LIFETIME  EXPERIMENT 

5.1  5.2  MeV  protons  from  the  University  of  Alberta 

van  de  Graaf  accelerator  were  used  to  bombard  the 

5  2  5  2 

Cr  target  and  initiate  the  Cr  (p,  ny)  Mn 
reaction.  In  order  to  obtain  a  reasonable  coinci¬ 
dence  rate  and  still  retain  good  energy  resolution  in 

the  Ge  (Li)  a  count  rate  in  the  Ge  (Li)  of  approximately 

3  -1 

8  x  10  sec.  was  maintained  throughout  the  experiment 
requiring  a  D.C.  beam  intensity  of  70  nA.  Coincident 
gamma  rays  were  observed  from  about  600  KeV  in  energy 
to  3100  KeV.  A  one  hour  singles  run  was  taken  at  an 
angle  of  90°  to  determine  the  zero-shifted  gamma  rays. 
The  Doppler  shifted  gamma  rays  were  obtained  with  a 
five  hour  coincidence  run  at  zero  degrees.  Using  an 
auxiliary  kicksorter,  the  SCA  windows  set  on  the  TAC 
spectra  were  periodically  checked  for  evidence  of  drift¬ 
ing  (see  figure  10  and  figure  11) . 

5.2  Results 

The  following  gamma  rays  we re  observed  in  the 
coincidence  and  the  singles  spectra: 


Figure  10 .  This  is  the  neutron  time-of- 
flight  spectrum  from  the  TAC  using  coin¬ 
cident  gamma  rays  as  the  zero  time  mark 
(large  peak) .  To  the  left  of  this  peak 
are  a  large  number  of  very  small  peaks 
forming  the  neutron  groups.  The  time 
resolution  is  not  good  enough  to  resolve 
a  particular  group  responsible  for  popu¬ 
lating  a  level  in  the  residual  nucleus. 
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Figure  11.  The  neutron  identifier  time 
spectrum  from  the  TAC  is  shown  and  results 
from  the  different  amounts  of  "slow"  com¬ 
ponent  in  gamma  ray  and  neutron  scintilla¬ 
tions  in  the  organic  liquid  detector. 

The  neutron  counts  have  been  raised  by  a 
factor  of  10  for  illustrative  purposes. 


LU 

cr 

3 

O 

Ll 


CHANNELS 


Figure  12.  The  portion  of  the  decay 
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scheme  of  Mn  considered  for  this 
experiment  is  shown  here. 
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1440  KeV 

(1440  - 

0) 

1288 

(1288  - 

0) 

910 

(1288  - 

— »  378) 

and  figure  12  shows  the  level  scheme  considered. 

Full  kinematic  shifts  were  determined  for  the 
three  gamma  ray  energies: 

910  KeV  3.003  KeV 

1288  4.250 

1440  4.694 

The  gain  shift  was  determined  by  calculating 
the  centroids  of  the  ^Co  source  peaks.  The  centroids 
obtained  for  the  ^Co  1.333  MeV  line  were: 

0°  centroid  (channel)  =  880.274  ±  0.08 
90°  centroid  =  878.736  ±  0.006  . 

The  difference  was  found  to  be  1.538  ±  .086  channels. 

6  D 

For  the  Co  1.172  MeV  line  the  centroids  were: 


j 
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0°  centroid  (channel)  =  774.943  ±  .08 
90°  centroid  =  773.410  ±  .006. 


Here  the  difference  was  1.533  ±  .086  channels.  Thus 
the  0°  spectrum  was  shifted  upwards  by  1.535  ±  .17 
channels.  Also  the  energy  calibration  for  both 
spectra  using  the  source  lines  was  obtained  by 
observing  that 


105.33  channels  =  160  KeV 

or  1  channel  =  1.517  KeV. 

The  same  centroid  analysis  was  applied  to  the 
three  gamma  ray  peaks. 

The  1288  — »  0  transition  (figure  13) 

0°  centroid  =  602.090  ±0.21 


600.468  ±  0.09 


Figure  13  and  Figure  14.  The  background 
subtracted  peaks  from  the  0°  and  90°  runs 
for  the  three  transitions  are  shown.  The 
background  for  the  0°  peaks  amounts  to  an 
almost  constant  10  counts  per  channel. 

For  the  1288  — *  0  and  1440  — *  0  90°  peaks, 

the  background  is  about  3000  counts  per 

channel  and  for  the  1288  - *  378  90°  peak 

it  is  6900  counts  per  channel.  The  un¬ 
corrected  shifts  are  also  indicated. 
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FIGURE  14 


Figure  15.  The  F(t)  function  is  plotted 

against  the  lifetime  x  (in  p.s.)  for  the 
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case  of  Mn  recoil  nuclei  slowing 
down  in  the  Cr  backing.  The  following 
experimental  points  with  F(x)  errors  are 
are  shown: 

(1)  1288  — >  378  and  1288  — 0 
transitions ; 

(2)  1440  — *■  0  transition. 
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shift  =  1.622  ±  0.30 

shift  with  gain  correction  =  1.622  -  1.535 

=  0.087  ±  0.47  channels. 

The  observed  shift  was  then 
0.132  ±  .71  KeV  . 


F(t) 


observed  shift 
kinematic  shift 


0.044  ±  0.236  . 


Having  used  the  F(x)  curve  given  in  figure  15  the 
lifetime  for  this  level  is 

1  c  +°° 

X  ~  1 . 3  J_3  }P  •  S  • 

The  large  upper  and  lower  bounds  on  the  lifetime  result 
from  the  F(t)  error  since  for  this  transition 

0  <  F (t)  <0.28. 

In  cases  where  the  data  falls  in  the  region  of  every 
small  F (t)  such  that  the  experimental  error  is  larger 
than  F(t)  one,  in  practice,  accepts  the  lower  bound  on 
the  lifetime.  That  is, 

t  >  0 . 2  p .  s  . 


The  1288  — »  0  transition  (figure  13) 
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0°  centroid  =  851.795  ±  0.30 

90°  =  850.144  ±  0.12 

difference  =  1.651  ±  0.42  channels 
corrected  shift  =  0.116  ±  0.590  channels 

=  0 . 176  ±  0 . 894  KeV 

F(t)  =  0.042  ±  0.211. 

Thus , 

T  =  1  5  +“ 

T  ^  * D  -1.3 

or 

t  >  0.2  p . s . 

The  1440  — ->  0  transition  (figure  14) 

0°  centroid  =  947.574  ±  0.39 

90°  =  946.019  ±  0.06 

difference  =  1.555  ±  0.43  channels 
corrected  shift  =  0.020  ±  0.62  channels 

=  0.030  ±  0.940  KeV 


F(t)  =  0.006  ±  0.190. 


Thus , 


T 


8 


+  00 


-7.7 


or 


t  >  0.3  p . s 
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Transition  t  exptl.(p.s.)  Weisskopf (p . s . ) 


1288 

3/2"  - 

378 

5/2" 

-f  CO 

1.5 

-1.3 

110 

0.042 

(E2) 

(Ml) 

1288 

0 

+  00 

1.5 

19 

(E2) 

3/2  - 

7/2 

-1.3 

1440 

0 

+  oo 

8.0 

11 

(E2) 

11/2 

— *  7/2 

-7.7 

Table  1 

These  results  along  with  the  single  particle  or 

Weisskopf  estimates  of  the  lifetimes  are  summarized  in 

table  1.  For  the  1288  — >  378  transition  it  seems  as 

though  E2  radiation  is  enhanced  and  Ml  radiation  is 

retarded.  The  other  two  transitions  are  pure  E2. 

Furthermore,  since  the  Weisskopf  estimate  for  El  radia- 

-4 

tion  is  so  low  (retarded  by  a  factor  of  10  from  E2) 
compared  to  the  experimental  limit  for  the  1288  — *  378 
transition  and  since  the  other  transitions  are  pure  E2, 
the  parity  for  the  first  three  excited  levels  are 


verified . 
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5.3  Throughout  the  experiment  the  problem  of  low  coin¬ 
cidence  rates  was  apparent  and  resulted  in  the  poor 
statistics  of  the  0°  spectrum.  This  condition  may  be 
improved  first  by  taking  care  to  see  that  the  neutron 
detector  discriminator  level  is  set  just  above  the  self¬ 
triggering  point  since  its  inherent  neutron  threshold  is 
approximately  800  KeV,  and  a  level  much  above  this 
(~  1.5  MeV)  will  result  in  a  loss  of  high  energy  coincidence 
gamma  rays.  Also,  the  bias  should  be  applied  several  hours 
before  run  time  to  avoid  immediate  drifting  of  the  discri¬ 
minator  level.  The  excessive  gamma  ray  count  rate  in  the 

1  " 

neutron  detector  indicates  that  the  /q  lead  shielding  is 
not  thick  enough  to  reduce  the  flux  of  low  energy  gamma 
rays  and  should  be  increased  (perhaps  to  h" ) . 

Secondly,  difficulty  was  encountered  in  trying  to 
get  all  the  signals  going  to  the  coincidence  units  to 
arrive  at  the  same  time  since  adjustment  of  an  SCA  window 
disturbed  the  timing.  This  problem  can  be  minimized  by 
strobing  the  SCA's  instead  of  the  TAC ' s  as  shown  in 
figure  7.  Also,  the  neutron  window  on  the  n-y  time  spectrum 
should  be  set  as  close  to  the  gamma  ray  peak  as  possible 
to  include  the  high  energy  neutrons  that  populate  the 
lower  levels  of  the  residual  neucleus. 

Finally,  modifications  should  be  made  to  the  target 
pot  assembly  including  better  beam  collimators  and  the 
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electrical  isolation  of  the  target  from  the  target  pot 
so  that  the  current  indicator  actually  registers  a  beam 
on  the  target  instead  of  on  the  entire  assembly.  The 
only  indicator  of  a  beam  on  a  target  during  the  experi¬ 
ment  was  the  gamma  ray  ratemeter. 

In  conclusion,  the  adoption  of  these  improvements 
will  serve  to  increase  the  n-y  coincidence  rate  (as  has 
been  shown  in  a  later  experiment  where  the  ratio  of  gamma 
ray  to  neutron  counts  in  the  neutron  detector  has  been 
decreased  from  15  to  1  to  3  to  1) ,  thus  making  the  present 
system  employing  coincidence  methods  very  useful  in 
measuring  Doppler  shifts  for  DSAM  lifetime  studies. 


■ 
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Appendix  1 

The  Angular  Distribution  Function  F(<j>,s) 

Suppose  a  charged  particle  is  moving  in  a  medium 
containing  N  atoms  per  unit  volume.  Its  position  is 
given  by  x  and  its  direction  of  motion  is  given  by  the 

/s. 

unit  vector  <j) .  The  scattering  cross-section  per  unit 
solid  angle  for  the  particle  in  the  material  is 
u(|$-cj),|)  (that  is,  the  cross-section  for  the  particle 
changing  its  direction  from  $  to  the  angle  between 

these  two  unit  vectors  being  the  scattering  angle  9). 

The  distribution  function  of  the  particle  will  be 
denoted  by  f(x,<j>,s),  where  s  is  the  arc  length  traversed 
by  the  moving  particle. 

Assume  now,  that  at  some  position  x,  arc  length 
s,  and  direction  $,  the  particle  suffers  a  collision 
which  changes  its  direction  to  $ 1 .  The  distribution 
function  changes  from  f(x,$,s)  to  f(x,$',s).  Conser¬ 
vation  of  particles  demands  that  the  change  in  f  due 
to  the  change  in  direction  of  the  particle  plus  the  sum 
of  the  changes  in  f  as  functions  of  x  and  s,  must  be 
equal  to  zero.  Or,  symbolically: 

- 

x,s,<j> 


0 


(1) 


- 
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The  divergence  of  f  with  respect  to  the  first 
two  variables  will  clearly  be 

3f  * 

“■  +  (J>  •  V  f  .  (2) 

The  change  of  f  with  $  will  be  given  by  the 

integral  of  f(x,$',s)  -  f(x,(J>/s)  weighted  by  the  cross- 

/\  /\ 

section  for  scattering  from  <j>  to  (j> '  over  all  final 
directions  $  1 :  Thus 

+  $•  Vf  =  Nf  [f  (x,  $  '  ,  s) -f  (x,  $  ,  s)  ]  a  (  |  $-$  '  |  )  d$* 
a  s 

(3) 

which  is  of  course  the  diffusion  equation  (1)  with  the 
last  term  transposed.  The  boundary  condition  of  this 
equation  will  be  that  f  be  unity  at  the  origin: 

f(x,$,0)  =  6  (x)  6($)  .  (4) 

$  =  0  means  that  the  particle  at  the  origin  is  assumed 
to  be  moving  in  the  incident  beam  direction  (z-axis) . 
f  is  now  expanded  in  terms  of  the  spherical  harmonics 
in  $: 

CO 

f(x,$,s)  =  I  f  (x,s)  Y  ($)  .  (5) 

£m 

=0 


Hence,  one  has  from  (3) 


I J 


?  .  ('#■’  je,' 
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3  f 


£m 


£m 


3  s  £m 


Ya mW  +  Z  (4>  Y£m^)).  Vf 


£m 


£m 


=  N  2  (^')-Y1)m($)]a(|?-$-|)d$'-  (6) 


£m 


£m  £m 


£m 


Multiplying  (6)  from  the  left  by  Y^_.  ($)  and 
then  integrating  with  respect  to  $  gives, 

Ts  +  $  V(?) 


=  N  [m  hm  SS  d<i> ' 


A  A 


=  N  Z  f  .  /  [P  (cos  0)  -  1]  6..  a(6)  d<)>  (7) 

Im  x'3  1 


using  the  addition  theorem  for  spherical  harmonics  and 


the  orthonormality  condition. 


Define 


A.  A  A  A 


Qi j  =  f  hjW  Ym($)  d$  . 


The  right  side  of  (7)  becomes 


2  TT  IT 

N  f .  .  /  /  [P •  (cos  0) 

30  0  0  1 


-  1]  a  (0)  sin0  d0  d0 


TT 


=  2  tt  N  f .  .  f  [P.  (cos  0)  -  1]  a(0)  sin0  d0 

ID  0  X 


=  -  f .  .  K. 
ID  i 


(8) 


. 


• 
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where  is  equal  to  -2  ttN  times  the  integral.  Thus 
(7)  becomes 


3f .  • 

— r-U-  +  K.  f.  .  = 
9  s  l  i  j 


-  Z  V  f 
£m 


£m 


Q 


£m 

ij 


(9) 


the  boundary  condition  of  which  is  given  by 


fio(^°)  =  fi(x,0)  =  6  (x)  Yio(0) 

=  «(x)  .  (10) 

j  (also  m)  has  been  set  equal  to  zero  since  the  cylin¬ 
drical  symmetry  of  the  problem  demands  this.  (Note 
that  because  m=j=0,  &=i) .  Thus,  one  has 


1 

/  $  P  ($)  d($)  =  0 

-1  36 


(11) 


since  the  integrand  is  an  odd  function.  Putting  these 
results  into  (9)  one  has 


9f  (x,s) 

-a—  +  K*  h(x's)  =  0  •  (12) 

Define 

* 

/  f£(x,s)  dx  =  F^(s)  • 

Integrating  (12)  over  all  space  will  then  give 

3*ys) 

a  s 


t  K,  FA(S) 


0 


(13) 


. 
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the  solution  of  which  is  given  by 


1"  s 

F£(s)  =  SXp  g  K£  ds) 


(14) 


where  F^(s)  satisfies  the  boundary  condition  (10) 
through  the  relation 


F£(0)  =  /  f^(x,0)  dx  =  (~ )  2  6 (x)  dx 


,21+l.k 
[  4tt  ; 


Returning  to  (5)  one  finds  that 


F£($,s)  H  f  f 0  (x,$,s)  dx 


Z  /f  (x,s)  dx  (cos  0) 

£  ,m=0 


00 


z  ( 2 it+1 )  P.  (cos  0 )  exp  (-  /  K  ds) 
4r  %  0  t 


(15) 
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Appendix  2 


Derivation  of  cos  $ 

Recall , 

i  r  +  a  . _  x 

cos  <f>  =  - - 5-7;  (i) 

(1  +  2ra  +  r  )  2 
where  a  =  cos  0. 

The  generating  function  for  Legendre  polynomials 
is  given  by 

2  -L  °°  n 

(1  -  2ra  +  r  )  2  =  E  P  (a)  r  .  (2) 

n=0  n 

From  the  properties  of  the  Legendre  polynomials 
one  has  that 


(-l)n  Pn(a) 


(3) 


and  the  recursion  relation 


( 2n  +  1)  a  P  (a)  =  (n  +  1)  pn+p(a)  +  n  pn_1  (a)  •  (4) 


Thus , 


(1  +  2ra  -i-  r2)  ^  =  Z  (-l)n  pn  (a)  r11 

n=0 


(5) 


and 


98 


cos  <J>  =  (r  +  a)  E  (-l)n  P  (a)  r11 

n=0  n 


=  £  (~l)n  P  (a)  rn+1  +  £  (-l)n  a  P  (a)  rn 

n=0  n=0  n 


=  -  Z  ("l)n  P„  ,  (a)  r11  +  a  +  £  (-l)n  a  P  (a)  rn. 

n=l  n  n=l 


(6) 


Since  n  is  a  dummy  index  it  has  been  replaced  by  n-1  in 
the  first  term  causing  the  summation  to  start  from  n=l. 
But 

X  (-l)n  a  P  (a)  rn  =  E  (-l)n  Pn+1  (a)  + 
n=l  n=l 


+  0  £  P  n  (a) ]  rn 
2n+l  n-1 


(7) 


so  that 


cos  <j)  =  a  +  E  (-1)  {- 

n=l 


n  r  n+1  _  ,  x  ,  r  n  n  T,  /  \  t  n 

(a)  +  U-  7T  -  1]  p_  n  (a)  }r 


2n+l  n+1 


2n+l 


n-1 

(8) 


using  (6)  and  (7).  Or, 


cos  <j>  =  P,  (cos  0)  +  E  (-l)n  - U+i-  x 
L  n=l 


2n+l 


n 


x  [p  .(cos  0)  -  P  ,  (cos  0) ]  r  .  (9) 

n+1  1 


n-1 


r  <  1 
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For  r  >  1  the  generating  function  is  given 


by 


(1  -  2ra  +  r2)  =  Z  P  (a)  (i)n+1 

n=0 


and  the  rest  of  the  derivation  proceeds  in  an  entirely 
analogous  manner  and  gives 


cos  <j>  =  1  +  Z  (~l)n  [P^Ll  (cos  0)-P„_1(cos  0)]  x 


n=l 


2n+l  n+1 


n-1 


,L  n+1 
x  (— ) 

r 


r 


>  1 


(10) 


